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Abstract 



Quenched QCD simulations on three volumes, 8^x, 12^ x and 16^ x 32 
and three couplings, (3 = 5.7, 5.85 and 6.0 using domain wall fermions pro- 
vide a consistent picture of quenched QCD. We demonstrate that the small 
induced effects of chiral symmetry breaking inherent in this formulation can 
be described by a residual mass {m^es) whose size decreases as the separation 
between the domain walls (L^) is increased. However, at stronger couplings 
much larger values of Lg are required to achieve a given physical value of rrires- 
For (5 = 6.0 and Lg = 16, we find mj-cs/ms = 0.033(3), while for /3 = 5.7, and 
Lg = 48, mres/fng = 0.074(5), where rUg is the strange quark mass. These 
values are significantly smaller than those obtained from a more naive de- 
termination in our earlier studies. Important effects of topological near zero 
modes which should afflict an accurate quenched calculation are easily visible 
in both the chiral condensate and the pion propagator. These effects can be 
controlled by working at an appropriately large volume. A non-linear behav- 
ior of in the limit of small quark mass suggests the presence of additional 
infrared subtlety in the quenched approximation. Good scaling is seen both 
in masses and in f^^ over our entire range, with inverse lattice spacing varying 
between 1 and 2 GeV. 



Typeset using REVT^ 



1 



I. INTRODUCTION 



Since spontaneous chiral symmetry breaking is a dominant property of the QCD vacuum 
and is responsible for much of the low energy physics seen in Nature, having a first principles 
formulation of lattice QCD which does not explicitly break chiral symmetry has been an 
important goal. Both Wilson and staggered fermions recover chiral symmetry in the contin- 
uum limit but with these techniques the chiral and continuum limits cannot be decoupled. 
For the QCD phase transition, which is dominantly a chiral symmetry restoring transition, 
a formulation that is free of violations of chiral symmetry due to lattice artifacts, should 
give a phase transition more closely approximating that of the continuum limit. For the 
measurement of matrix elements of operators in hadronic states, a formulation that respects 
chiral symmetry on the lattice substantially reduces operator mixing through renormaliza- 
tion. Lastly, since much of our analytic understanding of low-energy QCD is formulated in 
terms of low-energy effective field theories based on chiral symmetry, a lattice formulation 
preserving chiral sjTumetry allows controlled comparison with analytic expectations. 

Building on the work of Kaplan , who showed how to produce light chiral modes in a 
d dimensional theory as surface states in a + 1 dimensional theory, a number of attractive 
lattice formulations have been developed which achieve a decoupling of the continuum and 
chiral limits. Here we will use Kaplan's approach as was further developed by Narayanan 
and Neuberger and by Shamir 0. It is Shamir's approach, commonly known as the 
domain wall fermion formulation, which we adopt. (For reviews of this topic see Refs. [0-|Tl| 
and for more extensive recent references see Ref. |jl2[-) For a physical four- dimensional 
problem, the domain wall fermion Dirac operator, D, is a five-dimensional operator with 
free boundary conditions for the fermions in the new fifth dimension. The desired light, 
chiral fermions appear as states exponentially bound to the four-dimensional surfaces at the 
ends of the fifth dimension. The remaining modes for D are heavy and delocalized in the 
fifth dimension. 

An additional important feature of the domain wall fermion Dirac operator in the limit 

— > oo is the existence of an "index" , an integer that is invariant under small changes in 
the background gauge field. Here Lg is the extent of the lattice in the fifth dimension. This 
property, true for all but a set of gauge fields of measure zero, can be readily seen using 
the overlap formalism In the smooth background field limit, this index is the normal 

topological charge but, even for rough fields, it signals the presence of massless fermion 
mode(s) when non-zero. These zero modes can easily be recognized in numerical studies 
with semiclassical gauge field backgrounds [jl3|-|I^. 

These powerful theoretical developments in fermion formulations require additional study 
to demonstrate their merit for numerical work. For the case of domain wall fermions, a 
growing body of numerical results are available. Both quenched [p!^ , |l9| -|3l[| and dynami- 
cal pT|- p5| , P^ domain wall fermion simulations have been conducted and the domain wall 
approach is readily adapted to current algorithms for lattice QCD. (Much work is also be- 
ing done on the numerical implementation of the overlap formulation and its variations 
p6|-^.) A fundamental question, which is a major part of this paper, involves quantifying 
the residual chiral symmetry breaking effects of finite extent in the fifth dimension. 

Due to current limits on computer speed, some lattice QCD studies are only practical 
when the fermionic determinant is left out of the measure of the path integral. The resulting 
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quenched theory does not suppress gauge field configurations with hght fermionic modes, 
in contrast with the original theory where, for small quark mass, the determinant strongly 
damps such configurations. The measurement of observables involving fermion propagation 
through configurations with unsuppressed light fermionic modes can in principle lead to 
markedly different infrared behavior than that found in full QCD, in the limit of small 
quark masses. Domain wall fermions, which produce light chiral modes at finite lattice 
spacing and preserve the global symmetries of continuum QCD, should produce a well- 
defined chiral limit for full QCD. The central question addressed in this paper is whether 
a well-controlled chiral limit also exists within the quenched approximation. A thorough 
theoretical and numerical understanding of the quenched chiral limit is essential if the good 
chiral properties of domain wall fermions are to be exploited in quenched lattice simulations. 

Here we present results from extensive simulations of quenched QCD with domain wall 
fermions, primarily at two lattice spacings, ~ 1 and 2 GeV. Many different values for the 
fifth dimensional extent, Lg, and the bare quark mass, ruf have been used. Hadron masses, 

and the chiral condensate, (qq), are the primary hadronic observables we have studied. In 
calculating physical observables using domain wall fermions, four-dimensional quark fields 
q{x) are defined from the five- dimensional fields ^'(a;, s) by taking the left-handed fields 
from the four-dimensional hypersurface with smallest coordinate in the fifth dimension and 
the right-handed fields from the hypersurface with the largest value of this coordinate. We 
also present results from measuring the lowest eigenvectors and eigenvalues of the hermitian 
domain wall fermion operator. 

Here we list the major topics in each section of this paper. Section || defines our conven- 
tions and gives details of the hermitian domain wall fermion operator. Section ^ discusses 
our simulation parameters and fitting procedures and includes tables of run parameters and 
hadron masses for mj > 0.01. In Section ^ a precise understanding of how finite Lg effects 
enter (gg) is developed and measurements of (gg) which show the role of fermionic zero 
modes are reported. We study the pion mass in the chiral limit in Section |V|, which requires 
understanding zero mode effects. Section ^ contains two determinations of the residual 
chiral symmetry breaking for finite L^; one from measuring appropriate pion correlators and 
the other from the explicitly measured small eigenvalues and eigenvectors of the hermitian 
domain wall fermion operator. Our determination of /^r, an important check of the chiral 
properties of domain wall fermions, is discussed in Section |V11| , along with the scaling of 
hadron masses. 

Because of the length of this paper and the number of topics covered, we now give a 
brief summary of our major results, organized to correspond to the expanded discussion in 
Sections |rV||V][V3 and [VTI[ 

Zero mode effects in (gg): As already mentioned, the domain wall fermion operator D 
has an Atiyah-Singer index v for Lg ^ oo. However, in quenched QCD, v plays no role in 
the generation of gauge field configurations. For 7^ 0, both D and the hermitian domain 
wall fermion operator Dh |^ have zero modes. Since (gg) is an appropriately restricted 
trace of Dj^^ it should diverge as (|z/|)/mjV for small nif if the ensemble average of \h'\ is 
non-zero. Here V is the four-dimensional, space-time volume of the lattice being studied. 
For large but finite L^, the residual chiral symmetry breaking should cut off this divergence. 

Figure |l] shows (gg) versus the quark mass mj for ~ 1 GeV on two different volumes 
of linear dimensions of about 1.6 and 3.2 Fermi. A divergence for m/ — is clearly visible 
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on the smaller volume, but not on the larger. This is expected since /V should go as 
l/\/V and is clear evidence for unsuppressed zero modes in quenched QCD, first reported 
in Ref. [Q. Notice that there may be other problems with the chiral limit of (gg) that are 
masked by this divergence. 

The chiral limit of m^ : With this clear evidence for zero mode effects in (gg) , one might 
expect to see zero mode contributions in any quark propagator D^^{x, y) if at both x and y 
a single zero eigenvector has reasonable magnitude. For sufficiently large volume, needed to 
see asymptotic behavior in the limit of large \x — y\, there should be no zero mode effects. 
Our results for the zero mode effects on the pion mass are presented in Figure ^ which shows 

versus for 8^ x 32 lattices with Lg = 48 and Figure ^, where all the parameters are 
the same except that the volume was increased to 16^ x 32. The pion mass is determined 
from three different correlators which are each affected differently by zero modes. For the 
smaller volume, the pion masses measured disagree for small rrif, while they agree for the 
larger volume. 

Notice that on the larger volume shown in Figure |^, where zero-mode effects are not 
apparent, shows signs of curvature in mj with the three rrif = values lying below the 
extrapolation from larger masses. In addition, this simple large-mass linear extrapolation 
vanishes at a value of rrif that is more negative than the point mf + m^es = (shown in 
the graph by the star) also suggesting downward concavity. While the discrepancy between 
this X- intercept and the point rrif + nij-cs = 0, may be caused by O(a^) effects, we find a 
considerably larger discrepancy when making a similar comparison at /? = 6.0. Thus, we 
have evidence that does not depend linearly on m/ in the chiral limit. 

Determining the residual mass: In the limit of small lattice spacing, the dominant chiral 
symmetry breaking effect, due to the mixing between the domain walls, is the appearance of 
a residual mass, m^es in fhe low energy effective Lagrangian. The Ward-Takahashi identity 
for domain wall fermions |^ has an additional contribution representing this explicit chiral 
symmetry breaking due to finite Lg. Matrix elements of this additional term between low 
energy states determine the residual quark mass. Figure ^ shows our results for rrij-es for 
16^ X 32 lattices at /3 = 6.0 as a function of Lg. mres is clearly falling with Ls and reaches 
a value of ~ 2 MeV for Lg > 24. Our data does not resolve the precise behavior of mres 
for large Lg, but the very small value makes this less important for current simulations. 
A similar study on lattices with ~ 1 GeV but with larger Lg = 48 finds a value of 
rrires/ms = 0.074(5) or m^es ~ 8 MeV. 

We have also used the Rayleigh-Ritz method, implemented using the technique of Kalk- 



Reuter and Simma to determine the low-lying eigenvalues and eigenvectors for the 
hermitian domain wall fermion operator. The results exhibit the approximate behavior 
expected from low-energy excitations in the domain wall formulation. We use the resulting 
eigenvalues to provide an independent estimate of the residual mass which is nicely consistent 
with the more precise value determined from pseudoscalar correlators. 

Results for f^,, hadron masses and scaling: With this detailed understanding of the chi- 
ral limit of quenched lattice QCD with domain wall fermions, we have calculated /,r using 
both pseudoscalar and axial- vector correlators. The results for lattices with ~ 2 GeV 
are shown in Figure |, where good agreement between the two methods is seen. To do this 
comparison, the appropriate Z-factor for the local axial current must be determined and a 
consistent value for m^es must be known. The good agreement in the figure is a significant 



4 



test of these measurements as well as the chiral properties of domain wall fermions. We find 
very good scaling in the ratio jT^jrUp for ~ 1 to 2 GeV. For rnj^jmp scaling is within 
6%. We also find that — (gg) = (256(8)MeV)^ from our fi = 6.0 simulations. 



II. DOMAIN WALL FERMIONS 

In this section we first define our notation, including the domain wall fermion Dirac 
operator, and then derive the precise form of the Banks-Casher relation for domain wall 
fermions, to second order in the quark mass. In this paper, the variable x specifies the 
coordinates in the four- dimensional space-time volume, with extent L along each of the 
spatial directions and extent along the time direction, while s = 0, 1, . . . , — 1 is the 
coordinate of the fifth direction, with Ls assumed to be even. The space-time volume V is 
given hj V = L^Nf. The domain wall fermion operator acts on a five- dimensional fermion 
field, s), which has four spinor components. A generic four- dimensional fermion field, 
with four spin components, will be denoted by ip{x), while the specific four-dimensional 
fermion field defined from s) will be denoted by q{x). The space-time indices for 
vectors will be enclosed in parenthesis while for matrices they will be given as subscripts. 
Our general formalism follows that developed by Furman and Shamir 



A. Conventions 

The domain wall fermion operator is given by 

Dx,s;x',s' = ^s,s'DI^^, + dx^Dj-^g, (1) 



1 

+ (Ms - 4)5,,,. (2) 



D 



s,s' 



' (1 - 75)5s+l,s' + (1 + 75)^s-l,s' - 25s,s' 



2 
2 



(1 - 75)5s,L,-l5o,s' + (1 + 'l5)5s,ohs-l,s' ■ (3) 



Here, t/,,^ is the gauge field at site x in direction /i, and s and s' lie in the range < 
s,s' < Ls — 1. The five-dimensional mass, representing the height of the domain wall in 
Kaplan's original language, is given by M5, while mj directly couples the two domain walls 
at s = and s = Lg — 1. Since the light chiral modes should be exponentially bound to the 
domain walls, rrif mixes the two chiralities and is therefore the input bare quark mass. The 
value of M5 must be chosen to produce these light surface states and, in the free field case, 
< M5 < 2 produces a single fermion flavor with the left-hand chirality bound to s = and 
the right to s = L,, — 1. In order to use our pre-existing, high-performance Wilson fermion 
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operator computer program as part of our domain wall fermion operator, we have used the 
operator D above, which is the same as D^p of Ref. 0. 



Following Ref. we define the four-dimensional quark fields q{x) by 

g(x) = PL^(x,0) + PR*(x,L,-l) 
qix) = ^(x, - 1)Pl + ^(x, 0)Pr 



(4) 



where we have used the projection operators Pr,l = (1 =t 75)/2- Symmetry transformations 
of the five-dimensional fields yield a four- dimensional axial current 



Ls-l 

J2 sign 

s=0 



(5) 



Here 



^{x + fi, s){l + l,)Ul^,T^{x, s) - ^{x, s){l - j,)U,J-^{x + fi, s) (6) 



while the flavor matrices are normalized to obey Tr [f^t^) = 5°"^. The divergence of this 
current satisfies 



A^Alix) = 2mfmx) + 2Jl{x 



'5q\ 



(7) 



where A^f{x) = f\x) — f{x — /i) is a simple finite difference operator and the pseudoscalar 
density J^{x) is 



-^{x, - l)P^t-^{x, 0) + ^{x, 0)PRr^(a;, - 1) 
g(x)r75g(x). 



(8) 



This equation differs from the corresponding continuum expression by the presence of the 
J^q{x) term, which is built from point-split operators at Ls/2 and Ls/2 — 1 and is given by 



Jg^x) = -^(x, LJ2 - l)PLr*(x, LJ2) + ^(x, L,/2)PRr^(x, LJ2 - 1] 



(9) 



We will refer to this term as the "mid-point" contribution to the divergence of the axial 
current. 

This mid-point term adds an additional term to the axial Ward-Takahashi identities 
and modifies observables, like the pion mass, which are controlled by these identities. The 
Ward-Takahashi identity is 



A,{A;ix)Oiy)) = 2mf{r,{x)0{y)) + 2{J^x)0{y)) + i{5'^0{y)). 



(10) 



For operators, O made from the fields qiy) and it has been shown |^ that the Jg^ term 
in Eq. ^ vanishes for flavor non-singlet currents when — > oo. For the singlet current, this 
extra term generates the axial anomaly. The mid-point term represents the contribution of 
finite Lg effects on the low-energy physics of domain wall fermions. 
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B. Definition of the residual mass and the chiral limit 



For domain wall fermions, the axial transformation which leads to the Ward-Takahashi 
identity of Eq. ^ rotates the fermions in the two half-spaces along the fifth direction with 
opposite charges. For rrif = 0, the action is not invariant under this transformation due 
to the coupling of the left- and right-handed light surface states at the midpoint of the 
fifth dimension. This results in the additional term in the divergence of the axial current, 
as given in Eq. ^. In the — oo limit where the explicit mixing between the s ~ 
and s ~ — 1 states vanishes, this extra "mid-point" contribution will be zero and a 
continuum-like Ward-Takahashi identity will be realized. 

Since we must work at finite Lg it is useful to characterize the chiral symmetry breaking 
effects of mixing between the domain walls as precisely as possible. We do this by adopting 
the language of the Symanzik improvement program |^9|,0. Here we use an effective con- 
tinuum Lagrangian to reproduce to 0{a^) the amplitudes predicted by our lattice theory 
when evaluated at low momenta and finite lattice spacing. Clearly Cq is simply the con- 
tinuum QCD Lagrangian, while Ci will include the dimension-five, clover term: ipa^'^ipF^^^ 



51| . The chiral symmetry breaking effects of mixing between the domain walls will appear 
to lowest order in a as an additional, dimension three operator cx a~^e~"^'' ipip- This term 
represents the residual mass term that remains even after the explicit input chiral symmetry 
breaking parameter rrif has been set to zero. The next chiral symmetry breaking contribu- 
tion from domain wall mixing will be O(a^) smaller, appearing as a coefficient of order 
for the clover term. 

We define the chiral symmetry breaking parameter rrij-es so the complete coefficient of 
the mass term in Co is proportional to the simple sum ruf + m^es- While this is a precise 
definition of mj-es, valid for finite lattice spacing, a precise determination of mj-es in a lattice 
calculation will be impeded by the need to quantitatively account for the additional chiral 
symmetry breaking effects of terms of higher order in a. 

Close to the continuum limit, for long distance amplitudes, the Ward-Takahashi identity 
given in Eq. |lO|must agree with the corresponding identity in the effective continuum theory. 
Thus, for the non-singlet case, the sum of the first two terms on the right-hand side of Eq. ^ 
must be equivalent to an effective quark mass, meff = mj + m^esi times the pseudo-scalar 
density Jg . Thus, the residual mass, m-^es appears in the low energy identity: 

where this equality will hold up to O(a^) in low-momentum amplitudes. 

Thus, close to the continuum limit, mres in Eq. |Tl] is a universal measure of the chiral 
symmetry breaking effects of domain wall fermions for all low energy matrix elements, with 
corrections coming from terms of higher order in the lattice spacing. However, away from 
the continuum limit the O(a^) terms may be appreciable. In addition, if there are high 
energy scales entering an observable, such a low energy description is not valid and the 
explicit chiral symmetry breaking effects of finite Lg can be more complicated than a simple 
additive shift of the input quark mass by mres- 

Many aspects of the chiral behavior of the domain wall theory can be easily understood 
by reference to the more familiar Wilson fermion formulation. For finite Lg the domain wall 
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formulation can be viewed as an "on- and off-shell improved" version of Wilson fermions. The 
low energy effective Lagrangian for domain wall fermions is the same as that for the Wilson 
case except the coefficients of the chiral symmetry breaking terms are expected to decrease 
exponentially with Lg. Viewed in this way, one might expect to achieve a vanishing pion 
mass by fine-tuning to a critical value, irtfc in very much the same way as one fine-tunes 
K to for Wilson fermions. As the above discussions demonstrates, rrifc = mres + O(a^). 
Just as in the Wilson case, this limit can be interpreted as the approach to the critical 
surface of the Aoki phase [pl,p^,EO]. 



C. The hermitian domain wall fermion operator 



A hermitian operator Djj can be constructed from D through 



(12) 



where {R^)ss' = Ss,Ls-i-s' is the reflection in the fifth dimension around the five- dimensional 
midpoint, s = {Lg — l)/2. Writing out Dh gives 



Dh = 75 D^^^, 6s+s',Ls-i 

+ 75 [Pl Ss+s',Ls + Pr ^s+s',Ls-2 — Ss+s',Ls-l 



while as an explicit matrix in the s, s' indices: 



D 



H 



j,Pn 75(^11 - 1) \ 
75P,?75p"-l) 75i^L 



j,Pn 75(^11 - 1) I^Pl 
V75PII-I) 75PL 



-ruf-f^PR J 



The eigenfunctions and eigenvectors of Dh will be denoted by 

Dh^a„ = Ah^^^Ah 
with the five-dimensional propagator given by 

^'A^(X,S)^L 75 {R5)s,s' 



0(5) ^ 



A 



H 



(13) 



(14) 



(15) 



(16) 



(Grassmann variables in the Euclidean path integral will be denoted by \1/ and while the 
eigenfunctions of Dh will be denoted \l/t and Ku-) 



We will find it convenient to define three additional matrices 



{^b)s-^' = ^....9' sign ( ^\ ^ - s 



(17) 
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Qi5 = Pl Ss,o 5s', + Pr 5s,l,~i (18) 

and 

QS^^ = Pl 5s,Ls/2 Ss',Ls/2 + Pr ^s,l,/2-i <5s',l,/2-i- (19) 
The transformation which generates the current in Eq. ^ is 

^ ^^expi-ia^eV^). (20) 

The matrices Q^'^^ and Q^^^^ are the two parts of which correspond to terms in D = 
j^R^Dh which are not invariant under the transformation in Eq. |2y. The matrix Q^'"^ 
underhes the exphcit mass term and, in the original operator D, exphcitly mixes the s = 
and s = Ls — 1 walls. Likewise, the matrix Q^^'^^ is a "mid-point" matrix with non-zero 
elements only in the center of the fifth dimension. It represents the component of the 
operator D which connects the left and right half regions. These two contributions provide 
the terms on the right hand side of Eq. ^ and one easily finds 

{r5,D/f} = 2m/Q(-) + 2Q('"P). (21) 

Since it is expected that there are eigenvectors of Dh which are exponentially localized on 
the domain walls, we see that with rrif = and the limit ^ oo taken, Dh anticommutes 
with Fs in the subspace of these eigenvectors. This is the property expected for massless, 
four-dimensional fermions in the continuum in Euclidean space. 

Using the matrix, Q^'^\ we can write a simple form for the four- dimensional chiral con- 
densate, (qq) 

- im) = -T^J2i'^ix)q{x)) (22) 
^ ^(^(x,s)(i?5g^-)),,,,v[/(a;,s')) (23) 



12V 



X 



12^ Xx^X. 



^ {Ah\i,Q^-^\Ah) \ (25) 



12V A 



H 



where in the last line a bra/ket notation has been used. The large angle brackets indicate 
the average over an appropriate ensemble of gauge fields. 

We define the pion interpolating field as 7r"(x) = iq{x)t°''^^q{x) and then find that the 
pion two-point function is given by (no sum on a) 
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Note that the generators, do not appear in the spectral sum, since they merely serve to 
specify the contractions of the quark propagators and that 7r"(x) = iJ^{x). To investigate 
the extra term in the axial Ward-Takahashi identity, Eq. we will also have need to 
measure the correlation function between interpolating pion fields defined on the domain 
walls and the mid-point contribution to the divergence of the axial current, Jg^. We define 
a mid-point pion interpolating field by T^^^p^ix) = iJ^g{x) and the spectral decomposition 
for the correlator between interpolating pion operators on the wall and the midpoint is 

\Ah,A'^ ^^H^^H I 

We define a local axial current as A'^^{x) = q{x)t""-ff^'~f5q{x) and note that it is different 
from defined in Eq. ^ The two-point function of the zeroth component of this current, 
{Aq(x) Aq(0)) , has a form similar to Eq. with a factor of 70 multiplying each Q^'^^ and 
an overall minus sign. Finally, our scalar density is a{x) = q{x)q{x) and the connected 
correlator {a{x)a{0))c also has the form of Eq. ^ with a factor of 75 multiplying each Q*-™-* 
and an overall minus sign. 



III. HADRON MASSES FOR m/ > 0.01 

In this section we present the results for m^, rUp and ttt-at obtained for reasonably heavy 
input quark mass, m/ > 0.01 where the lower limit corresponds to mquark ~ ^strange/4. The 
more challenging study of m^r for m/ — > is described later, in Section |V|. This section is 
organized as follows. We begin by describing the Monte Carlo runs on which the results in 
this paper are based. Next the methods used to determine the hadron masses are discussed, 
both the propagator determinations and our fitting procedures. Finally, we present the 
results of those calculations for the easier, large mass case, rrif > 0.01. 



A. Simulation summary 

The results reported in this paper were obtained from ensembles of gauge field configu- 
ration generated from pure gauge simulations using the standard Wilson action |Q at three 
values of the coupling parameter, /5 = G/g"^: 5.7, 5.85 and 6.00. Thus, these ensembles follow 
the distribution, exp {6/(7^ X^-p trf/p} where the sum ranges over all elementary plaquettes 
V in the lattice and U-p is the ordered product of the four link matrices associated with 
the edges of the plaquette V. Some of the /? = 5.7 simulations and a portion of those at 



/3 = 5.85 were performed using the hybrid Monte Carlo algorithm |53]. These runs were 
performed on an 8^ x 32 space-time volume with a domain wall height M5 = 1.65. Each 
hybrid Monte Carlo trajectory consisted of 50 steps with a step size At = 0.02. These runs 
are summarized in Table |I[ In each case the first 2,000 hybrid Monte Carlo trajectories were 
discarded for thermalization before any measurements were made. After these thermaliza- 
tion trajectories, successive measurements of hadron masses and the chiral condensate, (qq) 
were made after each group of 200 trajectories. 
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A second set of simulations were performed using the heatbath method of Creutz [p4 



adapted for SU{3) using the two-subgroup technique of Cabibbo and Marinari and 
improved for a multi-processor machine by the algorithm of Kennedy and Pendleton ||56|| . 
The first 5,000 sweeps were discarded for thermalization. These runs are described in Table |I| 
where the values of M5 used are also given. Finally, the single j3 = 5.85 run with M5 = 1.9 
was performed using the MILC code ||5^. Here four over-relaxed heatbath sweeps |]58|j59[ 
with uj = 2 were followed by one Kennedy-Pendleton sweep, with 50,000 initial sweeps 
discarded for thermalization. 

A portion of the (3 = 5.7 masses described here appeared earlier in Ref. |2T] while the 
first of the (3 = 6.0 results appear in Refs. pO] and EH . 



B. Mass measurement techniques 

We follow the standard procedures for determining the hadron masses from a lattice 
calculation, extracting these masses from the exponential time decay of Euclidean-space, 
two-point correlation functions. In our calculation the source may take two forms. The first 
is a point source 

O'^ix) = q{x)erq{x) (28) 

which is usually introduced at the origin. The fiavor index a is introduced to make clear 
that we do not study the masses of fiavor singlet states. For the nucleon state we use a 
combination of three quark fields: 

Op(x) = eabcUaix)[ubix)C-f^dcix)] (29) 

where for simplicity we have written the source for a proton in terms of up and down quark 
fields, q = u and d. Here C is the 4x4 Dirac charge-conjugation matrix, e the anti-symmetric 
tensor in three dimensions and the color sum over the indices a, h and c is shown explicitly. 
Only these point sources are used in the /5 = 5.7 running. 

The second variety of source used in this work is a wall source. Such a source is obtained 



by a simple generalization of Eqs. ^ and ^ in which we replace the quark fields evaluated 
at the same space-time point x = (r, t) with distributed fields, each of which is summed over 
the entire spatial volume at a fixed time t. Gauge covariance is maintained by introducing 
a gauge field dependent color matrix V[U]{f,ti) which transforms the spatial links in the 
time slice t = ti into Coulomb gauge. Thus, to construct our wall sources we simply replace 
the quark field q{f,ti)c by the non-local field 

<l'"{U)c = E ^(^^ ti)c,c' q{r, U), (30) 

r 

where c and c' are color indices. We use these wall sources for the (3 = 5.85 calculations 
and a combination of both wall and point sources in the f3 = 6.0 studies. The use of wall 
sources for these weaker coupling runs is appropriate since the physical hadron states are 
larger in lattice units and better overlap is achieved with the states of interest by using these 
extended sources. 
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In all cases we use a zero-momentum-projected point sink for the second operator in the 



correlation function. This is obtained by simply summing the operators in Eqs. ^ and ^ 
over all spatial positions r in a fixed time plane t = tf. Thus, for example, we will extract 
the mass of the lightest meson with quantum numbers of the Dirac matrix F from the 
large tf — ti expression: 

r 

A similar equation is used for the nucleon correlation function except that the second expo- 
nent representing the state propagating through the antiperiodic boundary condition con- 
necting t = and t = Nf — 1 is reversed in sign and has exchanged upper and lower 
components for a spinor basis in which 7° is diagonal. 

For both the calculation of the quark propagators from which these hadron correlators are 
constructed and the evaluation of the chiral condensate, {qq) , we invert the five- dimensional 
domain wall fermion Dirac operator of Eq. |1], using the conjugate gradient method to solve 
an equation of the form Dy = h. This iterative method is run until a stopping condition is 
satisfied, which requires that the norm squared of the residual be a fixed, small fraction e 
of the norm squared of the source vector h. At the n*'* iteration, we determine the residual 
Vn as a cumulative approximation to the difference vector obtained by applying the Dirac 
operator to the present approximate solution ?/„ and h: = Di/n — h. We stop the process 
when |r„p/|/ip < e. 

For the calculation of (qq) we use e = 10~^ for the runs of Table H and e = 10^^ for 
those in Table |T|. For the computation of hadron masses in the runs of Table | we use 
e = 10^^ when Lg has the values 10, 16, 24 and 48, the condition e = 10^^'^ for the case 
Lg = 32. For the hadron masses computed in the runs in Table |I| we used e = 10^^ for 
f3 = 5.7 and 6.0, and e = 10"'' for f] = 5.85. Tests showed that zero-momentum projected 
hadronic propagators eight time slices from the source, calculated with a stopping condition 
of 10~^, differed by less than 1% from the same propagators calculated with a stopping 



condition of 10 for rrif > 0.01 For a quark mass rrif = 0.01 and a 16 x 32 volume 



with Ls = 16, typically ~ 1,500 conjugate gradient iterations were required to meet the 
stopping condition. For our very light quark masses {rrif < 0.001) up to 10,000 iterations 
were required for convergence. 

The final step in extracting the masses of the lowest-lying hadron states from the ex- 
ponential behavior of the correlation functions given in Eq. ^ is to perform a fit to this 
exponential form over a time range chosen so that this single-state description is accurate. 
Choosing tj = 0, we use the appropriate tf Nf — tf symmetry of Eq. ^ to fold the 
correlator data into one-half of the original time range < tf < Nf. We then perform a 



single-state fit of the form in Eq. |3T| for the time range tmin ^ tf < tmax < Nt/2. Typically 



^max is simply set to the largest possible value, tmax = 



^For the P = 5.7 runs we used smaller values of tmax for the vr and p fitting, typically 12 or 
14, in order to avoid the effects of rounding errors. These finite-precision errors, caused by a poor 
choice of initial solution vector, were seen at the largest time separations for the very rapidly falling 
propagators found at this strong coupling. 
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The lower limit, tmin, is decreased to include as large a time range as possible so as 
to extract the most accurate results. However, tmin must be sufficiently large that the 
asymptotic, single-state formula in Eq. ^ is a good description of the data in the time 
range studied. These issues are nicely represented by the effective mass, rriesit), with the 
parameters A and m = rricsit) in Eq. |3T| determined to exactly describe the hadron correlator 
at the times t and t + 1. To the extent that mcs{t) is independent of t, the data are in a 
time range which is consistent with the desired single state signal. As an illustration, this 
effective mass is plotted in Figure ^ for the vr, p and nucleon states in the 16'^ x 32, j3 = 6.0, 
Lg = 16, nif = 0.01 calculation. Good single-state ffis are easy to identify from the plateau 
regions for the case of and rrip. For the nucleon the rapidly increasing errors at larger 
time separations for this relatively light quark mass make it more difficult to determine a 
plateau. Better nucleon plateaus are seen for larger values of m/. 

The actual ffis are carried out by minimizing the correlated to determine the particle 
mass and propagation amplitude. We then choose tmin as small as possible consistent with 
two criteria. First, the ffi must remain sufficiently good that the P^r degree of freedom 
does not grow above 1-2. Second, we require that the mass values obtained agree with those 
determined from a larger value of tmin within their errors. 

In order to keep the ffiting procedure as simple and straight forward as possible, we 
choose values for t^[^ which can be used for as large a range of quark masses, domain wall 
separations and particle types as possible. Given the large number of Monte Carlo runs and 
variety of masses and Lg values it is possible to employ an essentially statistical technique 
to determine tmin- In choosing the appropriate tmm we examine two distributions. The ffist 
distribution is a simple histogram of values of x^/dof obtained for all quark masses and a 
particular physical quantum number. We require that for our choice of t^am, this distribution 
is sensibly peaked around the value 1 or lower. An example is shown in Figure |^ for the 
(3 = 6.00 p mass determined from a wall source for three values of tmin: 5, 7 and 9. 

In the second distribution we ffist determine a ffited mass mj(t) and the corresponding 
error cTj(t) for the state i, where the lower bound on the ffiting range is given by t. We then 
choose a t' > t and examine a measure of the degree to which mi{t) and mi{t') agree. The 
measure we choose is 

In Figure ^ we show the distribution of values of 6i{t', tmin) for the p meson for all t' > tmin 
and three choices for tmin: 5, 7 and 9. The distributions include p mesons with all values of 
rrif > 0.01 and all values for Lg used in the calculations. 

In our sample Figure 0, we have a reasonable distribution of x^/dof values for all three 
choices of t^i^ with only a slight improvement visible as tmin increases from 5 to 9. Likewise 
the distribution of mass values found at t' > t^i^ is in reasonable agreement for each value 
of train with a slight bias toward larger values being visible at the lowest value tmin = 5. 
Examining this figure and corresponding figures for the vr, for our quoted masses, we chose 
^min = 7 for these states. The fact that Figure |^ does not sharply discriminate between these 
three possible choices of tmin implies that we will get essentially equivalent results from each 
of these three values. 
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Our choices of tmin are as follows. For (3 = 5.7, where only point sources are used, tmin 
was chosen to be 7 for the tt, p and nucleon. For (3 = 5.85, hadron masses were determined 
only from the doubled 12^ x 64 configurations using wall sources and the value tmin = 6 for 
the p and 7 for the vr and nucleon. Finally for f] = 6.0 the most accurate mass values were 
determined using wall sources and it is these mass results which we quote below. Here tmin 
was chosen to be 7 for the vr and p and 8 for the nucleon. We were able to extract quite 
consistent results with larger errors using point sources. Here the needed value of t^in was 
10 for the vr and p and t^[^ = 8 for the nucleon. Finally, the errors are determined for each 
mass by a jackknife analysis performed on the resulting fitted mass. 



C. Hadron mass results 

The hadron masses that result from the fitting procedures described above are given in 
Tables [rV| - pCV| . Omitted from this tabulation are the masses for the more difficult cases 
ruf = 0.0 and 0.001 which are discussed later in Section 0. In each case the pion mass was 
determined from the {Aq{x)Aq{0)) correlator. While the results presented in these tables 
will be used in later sections of this paper, there are some important aspects of these results 
which will be discussed in this section. In particular, the dependence on volume and the ruf 
dependence of the p and nucleon masses will be examined. 

We begin by examining the dependence of the p and nucleon on the input quark mass. 



nif. In Figures ^, ^ and |T0| we plot the p and nucleon masses as a function of m/. As 
the figures show, each case is well described by a simple linear dependence on rrif. The 
data plotted in these figures appear in Tables |V11| pCll| and pClll| , respectively. Also plotted 



in Figure |T0] are our results for nip with non-degenerate quarks. The coincidence of these 
two results implies the familiar conclusion that to a good approximation the meson mass 
depends on the simple average of the quark masses of which it is composed. For simplicity in 
obtaining jackknife errors, we have included in these linear fits only that data associated with 
ensembles of configurations on which all relevant quark mass values were studied. Added 
configurations where only particular quark masses had been evaluated were not included. 

A simple linear fit provides a good approximation to all the masses considered in this 
section, in particular for mf > 0.01. In Table PCVI| we assemble the fit parameters for the 
/3 = 5.7, 8'^ X 32 masses, while Tables pCVll| , pCVlll| and pClX| contain the fit parameters for 



the p = 5.7, 16^ X 32, /3 = 5.85, 12^ x 32 and p = 6.0, 16^ x 32 calculations, respectively. 
The parameters presented in these three tables were obtained by minimizing a correlated 
which incorporated the effects of the correlation between hadron masses obtained with 
different valence quark masses, ruf, but determined on the same ensemble of quenched gauge 
configurations. The errors quoted follow from the jackknife method and the small values of 
X^/dof shown demonstrate how well these linear fits work. Because of the visible curvature 
in the pion mass for our f3 = 5.7 and 6.0 results, the linear fits for were made to the 
lowest three mass values. For the p and nucleon and all three masses at /3 = 5.85 we fit to 
the masses obtained for the full range of m/ values. 

Next we consider the effects of finite volume by comparing the 8'^ x 32 and 16^ x 32, 
volumes used in the P = 5.7, Lg = 48 calculation. The value of m.,, = 0.383(4) found at the 
lightest tTLf = 0.02 mass value for the 16^ x 32 implies a Compton wavelength of 2.6 in lattice 
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units. This lies between 1/4 and 1/3 of the hnear dimension of the smaller lattice, suggesting 
that we should not expect large finite volume effects. This is borne out by comparing the 
data in Tables |V111] and ^ where the two sets of masses agree within errors. 



This apparent volume independence within our errors can be nicely summarized by 
comparing the coefficients of the linear fits of the p and nucleon. Writing the two a and h co- 
efficients from the tables as a pair [a,b], we can compare the 16^ x 32 values from Table PCVII 



[0.775(18), 2.20(7)] and [1.03(4), 4.13(17)] for the p and nucleon with the corresponding num- 
bers for the 8^ x 32 numbers from Table ^VI|: [0.790(13), 2.18(5)] and [1.13(4), 3.79(15)]. For 
m^r the results on the two volumes agree to within the typical 1% statistical errors. However, 
for the case of the p and nucleon masses, finite volume effects may be visible on the two 
standard deviation or 1-2% level for the more accurate masses obtained for rrif > 0.06 

Since in lattice units the p mass decreases by about a factor of two as we change (3 from 
5.7 to 6.0, the 16^ spatial volume used at /? = 6.0 should be equivalent to the 8'^ volume just 
discussed at /3 = 5.7. Thus, we expect that the p and nucleon masses that we have found 
on this 16^ volume will differ from their large volume limits by an amount on the order of 
a few percent while the finite-volume pion masses may be accurate on the 0.5% level. 

IV. ZERO MODES AND THE CHIRAL CONDENSATE 

A. Banks-Casher formula for domain wall fermions 

In the previous section, our results for quark masses rrif > 0.01 were given, where 
the smallest values of m/ gave rriT^/mp ~ 0.4. Since the domain wall fermion operator 
with ruf = should give exact fermionic zero modes as Lg oo, observables determined 
from quark propagators at finite Lg, when small quark masses are used, should show the 
effects of topological near-zero modes. For quenched simulations, where zero modes are not 
suppressed by the fermion determinant, these modes can be expected to produce pronounced 
effects. One important practical question is the size of the quark mass where the effects are 
measurable. To begin to investigate this we now turn to the simplest observable where they 
can occur, {qq). 

Before considering the domain wall fermion operator, we review the spectral decomposi- 
tion of the continuum four-dimensional, anti-hermitian Euclidean Dirac operator The 
eigenfunctions and corresponding eigenvalues of such an anti-hermitian operator satisfy 

(^(4) + ^)^^ = + ^)^^. (33) 

with A real and 



The naive lattice fermion operator^ = 'yu{Ux,fiSx+fi,x' — uj., Sx-ji^x') and the lattice staggered 



fermion operator have eigenvalues and eigenvectors which also obey Eq. 
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(We use A to label eigenf unctions and eigenvalues of the anti-hermitian operator, saving Xh 
for the "hermitian" case defined below.) In the continuum, the presence of zero modes is 
guaranteed by the Atiyah-Singer index theorem for a gluonic field background with non-zero 
winding number []62| , |63| . 

The four-dimensional quark propagator, S^^^, can be written as 



^(4) ^ Mx)i^iiy) 



""'^ "j" tX + m ^ ^ 

leading directly to the Banks-Casher relation (with our normalization for the chiral 
condensate) 

(37) 

where u is the winding number and p(A) is the average density of eigenvalues. For quenched 
QCD, p(A) has no dependence on the quark mass. For both quenched and full QCD, one 
expects that |z/| ~ y/V, as is the case for a dilute instanton gas model. Thus, zero modes 
lead to a divergent 1/m term in {qq) whose coefficient decreases as 1/y/V. (This contrasts 



with the behavior seen |31[] above the deconfinement transition where it can be shown that 
the 1/m term remains non-zero for quenched QCD in the infinite volume limit [Q.) Before 
discussing the results of our simulations, we first address how this simple expectation of a 
1/m term in (qq) due to zero modes should appear for the domain wall fermion operator. 

We will find it useful to compare the spectrum and properties of the hermitian domain 
wall fermion operator with the hermitian four- dimensional operator, Dff\ defined by 

Dg)=75(^(')+m). (38) 

The eigenvalues, Xh, and eigenvectors, ipx^, for this operator can be given in terms of 
A and ipx given above. If A = we immediately get an eigenvalue Xh = for the 
hermitian operator, and an eigenvector with the definite chirality -|-1 or —1. For A 7^ 0, the 
eigenvectors of D^^ are linear combinations of ipx and ip^x and the corresponding eigenvalues 
are A^ = ±VA^ -|- m^. Since (^^^^ -|- m)"^ = (D^-')"^75, we have 



{^ij) = j^{Tr{p' + m)-') (39) 



1 



12V 



Tr(75(/^2^)-^)) (40) 



1 / ^A„(a^)W5^A„(a^) \ (^^) 



12V 



Since for D^^ 



E<(^)75^Aja:) = {^ (42) 
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Eq. ^ also reduces to the Banks-Casher relation, Eq. |3^. For finite mass, the zero-mode 
hermitian eigenfunctions are chiral, while other eigenf unctions have a chirality proportional 
to the mass. This will be important in our comparisons with domain wall fermions. 

For large Lg, it is expected that the spectrum of light eigenvalues of the hermitian domain 
wall fermion operator, Dh, should reproduce the features of the operator -D^\ Since Dh 
depends continuously on rrif, for small rrif its ith eigenvalue must have the form 

= < + b'i^f + Cim} + ■■■. (43) 



To make a connection with the normal continuum form for the eigenvalues we reparameterize 

A% i as 

Al, = nliX^ + imj + 6m,f ) + .... (44) 

Here n^^i is an overall normalization factor and we have defined 6mi, which enters as a 
contribution to the total quark mass for the ith eigenvalue. For 6mi = —rrif, A^j is at its 
minimum. Modes which become precise zero modes when Lg oo will have non-zero values 
for Aj and 6mi for finite L^. We will refer to such modes as topological near- zero modes. 
From perturbation theory in ruf, one can easily see that 



dA 



HA 



dnif 



(A//,|75Q^"^|Ah,) (45) 



while the chain rule applied to Eq. ^ gives 

dAH,i{mf) ^ nlijmf + Sruj) 
dm f A H,i 

Combining this with Eq. |25| gives 



(46) 



which agrees with the Banks-Casher form, Eq. ^ with the addition of the i dependent 
mass contribution Srrii. Thus, the parameter Aj in Eq. ^ should be identified with the 
eigenvalues of the continuum anti-hermitian operator As indicated by Eqs. ^ and 



6mi should represent a contribution to the eigenvalue from the chiral symmetry breaking 
effects of coupling of the domain walls, present for finite Lg. 

These arguments show that the domain wall fermion chiral condensate will grow as l/frif 
for gauge field configurations with topology, provided Lg is large enough to make 6mi and 
Aj small. The continuum expectation of a l/mj divergence is modified at small by the 
non-zero values of 6mi and Aj for topological near-zero modes. For a single configuration, 
the precise departure from a l/rrif divergence is dominated by the eigenvalues with the 
smallest values for 6mi and A^; for an ensemble average, the departure from l/rrif behavior 
depends on the distribution of values of 6mi. With this understanding of (qq) for domain 
wall fermions, we turn to our simulation results. 
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B. Quenched measurements of (qq) 



In this section we discuss our results for (qq) for quenched QCD simulations with domain 
wall fermions. Tables |, || and |Tl| give details about the runs where (qq) was measured. 



The most important aspect of the run parameters is the small values for nif used, including 
TTLf = 0.0 where finite Lg keeps non-zero, allowing the conjugate gradient inverter to be 
used. Of course the number of conjugate gradient iterations becomes quite large. 

Equation ^ shows that we should expect large values for — {qq) for small m / for con- 
figurations with topological near- zero modes. Figure |lT| shows —{qq) for 8'^ x 32 lattices at 
P = 5.7 with both Lg = 32 and 48. The quark masses used cover the ranges — 0.04 and 
0.00025 — 0.008, defined in Table |T|. Both values for Lg show an increase in —{qq) for very 
small quark mass, an effect expected from the presence of a no n- zero value for (This 
effect was first reported for domain wall fermions based on quenched simulations done on 
8^ X 32 lattices with jS = 5.85, Mf = 1.65 and = 32 and hsted in Table | [p2l.) 



Motivated by the form of Eq. we have fit — {qq) to the following phenomenological 
form 

~ = ^ |"a1 + + ""l""^ ('^^^ 

where a_i, Oq, ai and dnii^-qq^^ are parameters to be determined, dmi^qq^ represents a weighted 
average of dnii over the eigenvalues which dominate (gg) for small mf. The measurements 
of (gg) for different values of m/ are strongly correlated, being done on the same gauge 
field configurations with, generally, the same random noise estimator used to determine (gg) 
for all the masses. The common noise source makes the signal for the 1/mj divergence 
particularly clean, since the overlap of the topological near- zero mode eigenvectors with the 
random source does not fluctuate on a single configuration. This strong correlation precludes 
doing a correlated fit of (gg) to m/, since the correlation matrix is too singular. Thus, the 
fits in this section are uncorrected fits of (gg) to m/. 



Table PCX| gives the results for fits to the form of Eq. ^ for our (3 = 5.7, 5.85 and 
6.0 simulations. All the fits have a value x^/dof less than 0.1, a consequence of doing 
uncorrelated fits to such correlated data. In Figure |ll], one sees that the fit represents 
the data quite well. Continuing with 8'^ x 32 lattices at (3 = 5.7, Table pCX| shows the 
fit parameters are very similar for Ls = 32 and 48, except for 5m^gq), which drops from 
0.0040(4) to 0.0017(2). This indicates a decrease in Srui as Lg increases. 

Figure |12| is a similar plot of —{qq) for 16^ x 32 lattices with P = 6.0 for Lg = 16 and 
24. The rise in (gg) for small rrif exhibits the same general structure as for the (3 = 5.7 
data in Figure |lT|, but the effect is larger. Here dmi^qq^ falls from 0.00056(3) for Lg = 16 to 
0.00011(1) for Lg = 24. 

To further demonstrate that the divergence for small is due to eigenfunctions of 



Dh that represent zero modes of a definite chirality. Figure ^ shows the evolution of both 
— (gg) (solid lines) and — (g75g) (dotted lines). These evolutions are for 16^ x 32 lattices at 
P = 6.0 with Lg = 16. Eigenfunctions with a positive chirality contribute equally to (gg) and 
(g75g), while negative chirality eigenfunctions contribute with an opposite sign to (g75g). 
The topological near-zero modes should be approximately chiral and, for smaller values of 
rrif, one see large fluctuations in (gg) and (g75g). Some of the fluctuations have the same 
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sign and some are of opposite sign. Thus, we have configurations with eigenfunctions which 
are very good approximations to the exact zero modes expected as — oo. 

As mentioned earher, /V should decrease with volume, with the asymptotic depen- 
dence given by \/\/V. To investigate this numerically, we have measured (gg) on both 
8^ X 32 and 16^ x 32 lattices at /3 = 5.7 and 5.85 with = 32 and show the /3 = 5.7 results 
in Figure j^. The graph clearly shows that the l/rrif divergence is drastically suppressed by 
the larger volume. The coefficient of the 1/m/ term falls from 6.0(6) x 10~^ to 2.5(4) x 10~^ 
as the volume is changed by a factor of 8. This may be somewhat misleading, since ^m^g^^ 
also changes by a factor of about 2, likely due to the phenomeno logical nature of the fit 
and the small effects of the l/nif pole for the larger volume. Putting aside this systematic 
difficulty, the 1/mj coefficient decreases by a factor of 1/a/5.8, showing the general behavior 
expected but not in precise agreement with the expected asymptotic form. For (3 = 5.85, 
where the physical size of the lattices is smaller, the 1/mf coefficient falls from 3.8(3) x 10~^ 
to 0.60(9) X 10~^, a factor of 6.3. We have not seen the expected l/y/V dependence for the 
l/fTLf coefficient, but it does decrease with volume in accordance with general ideas. It is 
possible that on the larger 16^ x 32 volume, the l/rrif rise is not large enough to allow its 
coefficient to be determined without systematic errors. 

Thus, we have clear evidence for topological near- zero modes in our quenched simulations 
using domain wall fermions. They are revealed through a large l/rrif rise in our values for 
— (qq) , the presence of configurations where (qq) and {q'y^q) are large and of opposite sign and 
the volume dependence of the coefficient of the 1/mj term. We have extracted a quantity, 
6m(^qq), from a phenomenological fit to (qq), which represents the effects of finite Lg on the 
eigenmodes with small eigenvalues which dominate (qq) for — >■ 0. Physical values for 
(qq) in the chiral limit, without the contribution of the topological near- zero modes, will be 
presented in Section |V11| . We now turn to a discussion of how these zero modes, and the 
expected light modes responsible for chiral symmetry breaking, are evident in measurements 
of the pion mass. 



V. THE PION MASS IN THE CHIRAL LIMIT 

For domain wall fermions with Lg = oo, the chiral limit is achieved by taking rrif = 0. 
For our quenched simulations at finite L^, we must investigate the chiral limit in detail to 
demonstrate that the changes from the — oo limit are under control and of a known 
size. As is discussed in Sec. [11 B|, for low energy QCD physics the dominant effect of finite 



Lg should be the appearance of an additional chiral symmetry breaking term in the effective 
Lagrangian describing QCD. This term has the form rrires q{x)q{x) and in the continuum limit 
its presence will make vanish at = — mres up to terms of order a^. Our investigation 
of the chiral limit is made more difficult since there are other issues affecting this limit, 
beyond having Lg finite. For domain wall fermion quenched simulations, the chiral limit 
may be distorted by: 

1. Order effects. Since we are working at finite lattice spacing chiral symmetry will not 
be precisely restored even for m/ + mres = 0. In particular, additional chiral symmetry 
breaking will come from the effects of higher dimension operators suppressed by factors 
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of a" for n > 2. Thus, we cannot not expect m^r to vanish precisely at the point 
TTLf = — mj-es, but perhaps at a nearby point, removed from — mres by a terms of O(a^). 

2. Finite Ls- The residual mass, mres, should represent the finite Ls effects for physics 
describable by a low-energy effective Lagrangian. However, there will be additional 
effects of finite Lg for observables sensitive to ultraviolet phenomena. Further, a quan- 
tity with sufficiently severe infrared singularity may show unphysical sensitivity to 
those Ls-dependent eigenfunctions \AH,i) (and the parameters n^^i, Xi and Srui of the 
previous section) with small eigenvalues An^i- 

3. Topological near- zero modes. The previous section has shown these dominate (qq) 
for small quark masses { rrif < 0.01 ) for the volumes we are using. From the 
Ward-Takahashi identity, these effects must also by present in the pion correlator 
(7r"(x)7r'^(0)). 

4. Finite volume. For staggered fermions, where the remnant chiral symmetry at finite 
lattice spacing requires = when the input quark mass is zero, the finite volumes 
used in simulations have been seen to make non-zero when extrapolated to the 
chiral limit from above [|66| , |67| . Such an effect may also be expected to occur for 



domain wall fermions. 

5. Analytic results argue for the presence of "quenched chiral logs" with the dependence 
of ml on the quark mass in quenched QCD different from that of full QCD p8|-[70[]. 



In this section we study the pion mass in the limit of small quark mass. Demonstrating 
consistent chiral behavior for the pion mass in the limit mf+irires — is a critical component 
in establishing the ability of the domain wall fermion formalism to adequately describe chiral 
physics. If we discover that the limit ruf + m^es — > is obscured by large 0{a?) effects or 
large violations of chiral symmetry caused by unanticipated propagation between the domain 
walls, little may be gained from this new formalism. For the p and nucleon masses reported 
in Section |T|, the masses were shown to be well fit by a linear dependence on the input quark 
mass, mf. Any possible non-linearities are not resolvable within our statistics. For the pion, 
the statistical errors for these values of mj are smaller and we have also run simulations at 
smaller values for mf so we might hope to learn more about this important quantity. We 
begin by investigating the effects of topological near-zero modes on the pion. 



A. Topological near-zero mode effects on the pion: analytic considerations 

We have seen that topological near-zero modes dominate (gg) for small mf and, by con- 
tinuity, they will also alter the value for (gg) determined with larger quark masses. Through 
the Ward-Takahashi identity, these modes also appear in the pion correlator, (7r"(x)7r'*(0)) 
and therefore can enter in the determination of the pion mass in a lattice simulation. Alter- 
natively, the axial- vector correlator can be used to measure the pion mass and the zero modes 
may affect this correlator differently. It is vital to understand the role of these topological 
near-zero modes, since a study of the chiral limit of depends on an accurate measurement 
of the mass of the pion state. In this section we will study the way in which topological 
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zero-modes might be expected to effect pion correlation functions for the continuum theory 
using our results as a guide to the study of the domain wall amplitudes. 

Before proceeding, we first establish our notation for susceptibilities and the integrated 
Ward-Takahashi identity. In general we define 

XCD ^ ^E{C{x)Dm (49) 
where C and D are any two hadronic interpolating fields. In particular 

x.. = ^EK(^K(o)) (50) 

X^A = ^E(^S(^)^S(0)) (51) 

X.(.p). = ^EKmp)(^)^"(0)) (52) 

^'^ X 

where no sum over a is intended and the factor of 1/12 has been introduced to maintain 
consistency with our somewhat unconventional normalization for the chiral condensate given 
in Eq. ^ Then the Ward-Takahashi identity, Eq. |10|, with 0(0) = 7r"(0) and summed over 



X becomes 



mfXnTT + X7r(^p)7r = (53) 



which we will refer to as the integrated Ward-Takahashi identity. 

We first consider Eq. for large Lg, where we should recover the continuum version 



of the identity. To simplify the presentation, we start with the notation of Section [IV A 
for the continuum four-dimensional anti-hermitian Dirac operator. We immediately deduce 
from Eq. ^ that a 1/m divergence in (qq) from topological zero modes dictates a 1/m^ 
divergence in x-Kn- In addition, Xmr should have a 1/m divergence for large volumes from 
the pion pole and, as we will see below, there can also be a 1/m pole from topological zero 
modes. However, the volume dependence of these various pole terms should be different. 
Pole terms from topological near- zero modes should have a coefficient which is OiV^^^"^) for 
large volumes, while the 1/m term from the pion pole should be volume independent. 
Thus, we expect 

X.. = T7^^ + — + 0(m°). (54) 



m^ m 



The coefficients a_2 and a_i should become volume independent in the infinite volume limit. 
However, the "pion pole" piece, a_i, may contain an additional term arising from 

zero modes. Note, a particular order of limits must be understood when interpreting Eq. |5^. 
One expects that the usual relation oc m will hold only when m ^ 1/^1(^9)1 [O- 



Although this prevents our taking the m — limit of Eq. 5|, it is fully consistent with 
the domain m oc 1/V^/^ where the 1/m^ term in Eq. ^ may be as large as or much larger 
than the conventional 1/m term coming from the pion. For domain wall fermions at finite 
Lg these pole terms will be rendered less singular by the presence of the 6mi terms in the 
eigenvalues for Dh- 
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Lattice measurements of the pion mass come from the exponential decay of a correlator 
like (7r"(a;)7r°(0)) in the limit of large |x|. Having examined the zero mode effect in the 
somewhat simpler susceptibilities, we will now investigate the topological zero mode con- 
tributions to two-point functions from their spectral decomposition to understand how zero 
modes can distort measurements of the pion mass. We have 

(7r'^(x)7r'^(0)) = (Tr[ 5^75 5^5 75]) (55) 

^ / ^l(a;)^v(x) 4(0)^a(0) \ . . 

(-^A + m)(zA' + m) /' ^ ^ 

First we consider the terms in the sum where both A and A' are zero. This gives a 1/m? 
pole in (7r"(2;)7r'*(0)), provided the eigenfunctions in the numerator are non-zero at x and 0. 
(Integrating (7r"(x)7r"(0)) over x shows that these topological near-zero modes give the 1/m? 
contribution to x-k-k-) The terms in the sum where neither A or A' are zero should include 
the small eigenvalues which are responsible for the Goldstone nature of the pion. For large 
the total contribution to (7r"(a;)7r"(0)) from these modes should be proportional to 

|(0|7r(0)|7r)p r ^e-'"-!"!. (57) 



(Integrating over x gives another factor of m^r in the denominator, which produces the 1/m 
pion pole in Xtttt-) Lastly, the terms with either A or A' zero, but not both, can be written 
as 




A^ + m^ / X^^X^ + m?/ m ' 

where (qq) nz is the chiral condensate measured without zero mode contributions. Here we 
have used the symmetries in Eq. |3^ to combine the ±A terms in the sum over A and remove 
the term odd in iX. Since {qq)nz should be non-zero as m — *• 0, we see that the contribution 
to the correlator from terms with zero modes in one of the propagators can produce at most 
a 1/m pole term in (7r"(x)7r"(0)). 
Thus, we expect 

(..,.).»,o)> ^ -L (2^ + 2M^ , ^ . . . 



for small m. The first two terms represent the possible zero mode contributions. It is 
important to note that c_2(x,0) gets contributions from the modulus squared of the zero- 
mode eigenfunctions at the points and x, while c_i(x, 0) does not. In particular, for a 
configuration with a single zero mode, c_2(x, 0) is positive definite, being given by 

c.2{x,0) = V'/'{\M^)\' l^o(O)P). (60) 

Thus, one could expect c_2(a;, 0) to be a number of order the inverse of the mean zero- 
mode size squared, while c_i(a;, 0) could be much smaller due to the terms of differing sign 
appearing in the sum over eigenmodes. 
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For large enough only the true pion state should contribute. Such large |x| requires 
a correspondingly large V with the necessary suppression of zero modes. However, at a 
fixed separation in a finite volume and for simulations with small enough m, the physical 
pion contribution to (7r"(x)7r"(0)) can be completely negligible. For finite Lg, the domain 
wall fermion spectral form, Eq. |2^, gives the precise role of the topological near-zero modes. 
The double sum over A^^ and A'^j decomposes as we have done above for D^^^ and the 
dominant contribution of the topological near-zero modes enters as 1/(A^ -|- (m/ + 5mj)^) 
provided (A/^IQ^'^-'IA^) is well approximated by 0(1)5ah,a'j^- Thus, for Aj and Sirii small, 
there should be a region in rrif where (7r"(x)7r"(0)) displays a l/mj character. 

The pion mass can also be measured from the axial current correlator, {Aq{x)Aq{0)). The 
susceptibility for this correlator, xaa, is not constrained by the integrated Ward-Takahashi 
identity as is Xmr- However, there must be a pion pole contribution in addition to any zero 
mode terms. Therefore 

XAA ~ ^Ql^oWI^^I' + 0(m-") zero mode poles (61) 
mi 

TT 

= 0{m^) + 0{m~''^) zero mode poles (62) 

where we have used (0|y4o(0)|7r) ~ m^. The physical pion contribution is independent of m 
for small m, which is to be compared with the 1/m contribution of the physical pion to Xnn- 
We now turn to the question of the zero mode contribution. 

Once again we consider the Lg = oo case and use the notation for D^^\ For the axial 
vector correlation function, the spectral decomposition is 

(AS(x)Ag(O)) = - (Tr [ si'), 70 75 ^£ 7o 75 ] ) (63) 

^ ^i{x)io^x'{x) ^1,(0)7o^a(0) \ .g^. 

{-iX + m){tX' + m) /' ^ ^ 

The terms in the sum where both A and A' are zero modes vanish here, since the zero 
modes have a definite chirality and 70 couples different chiral components. (On a given 
configuration, all the exact zero modes must have the same chirality since exact zero modes 
can only occur through the index theorem.) Thus, there are no terms in {Aq{x) A'^{0)) . 
Note that a 1/m contribution to {Aq{x)Aq{0)) can appear from terms in the sum with either 
A or A' a zero mode (as we saw for (7r"(x)7r"(0))). The size of such a contribution depends 
on the matrix element of 70 between eigenfunctions. 

The terms with neither A nor A' a zero mode give the physical pion contribution, which 
should have the form 

ml ~ v^e-'"-!"!. (65) 



Thus, we expect that 



(A«(a:)A«(0)) = + di/2m'/' e"'"-!^! + • • • (66) 



d-i{x, 0) 

with other possible subleading terms from topological zero modes. As for the coefficient 
c_i(x, 0) in (7r'^(a;)7r"(0)), the coefficient d-i{x,0) above involves matrix elements between 
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different eigenf unctions and could be quite small from cancellations. Thus, even though 
in {Aq{x)Aq(0)) , the physical pion contribution can still be 0(m^/^) smaller than the zero 
mode contribution, the effects of zero modes in this correlator are likely suppressed by the 
smaller coefficient (i_i(a;, 0). 

To finish our discussion of the topological zero modes in correlators, we now examine 
the spectral form of {a{x)a{0))c, where the c subscript means that we only consider the 
connected part of the correlator. We find 

\t^' (-iX + m)(iX' + m) / ^ ' 

Since zero modes are eigenfunctions of 75, their contribution to the 1/m^ and 1/m terms in 
(7r''(a;)7r"(0)) and — (cr(a;)cr(0))c are equal. Thus, we have 

- (.(.).(0,). ^ ^ + ^-^) ^ e--N + . . . (68) 

for small m. By considering (7r"(x)7r"(0)) + ((t(x)(t(0))c, we obtain a two point function with 
no zero mode effects, but which contains both the physical pion and a heavier state from 
{a{x)a{{)))c- Thus, to reduce the effects of topological near-zero modes in this way requires 
that one works with correlators where the heavy mass ac state is present. 

To summarize this section, we have seen how the topological near-zero modes for domain 
wall fermions should enter the correlators which are used to determine the pion mass. For 
(7r''(a;)7r"(0)), there must be a l/roj contribution from near-zero modes, compared with 
the 1/ ^uif contribution expected from the physical pion. For {Aq{x) Aq{{))) ^ the topological 
near-zero modes can contribute a term of order l/ruf, while the physical pion should produce 
a ^ymf contribution. However, the coefficient of the term can be small. We also have 

pointed out that the volume dependence of the contribution of the topological near-zero 
modes to the correlator is different from the contribution due to the modes responsible 
for chiral symmetry breaking in QCD so that the zero-mode effects should vanish as the 
space-time volume increases. 

The above discussion explicitly addresses the behavior to be found in a chiral theory. 
Thus, it will apply to the domain wall theory in the limit — > 00. We might expect 
two sorts of modified behavior for a theory with finite L^. First, the chiral properties of 
the exact zero modes which eliminate the most singular terms from the {Aq{x)Aq{0)) and 
(7r"(x)7r"(0)) + {(t{x)(t{0))c will no longer be exact for finite Lg allowing more singular terms 
suppressed exponentially in to appear. Second the zero-mode singularities themselves 
may be softened by additional mass contributions to the denominators. We now turn to the 
results of our simulations. 



B. Topological near-zero mode effects on the pion: numerical results 

The first detailed studies of as rrif 0, done on 8^ x 32 lattices with (3 = 5.7 and 
a variety of values of L^, showed that ml{mf = 0) was not decreasing exponentially to 
zero as Lg — > 00, but rather seemed to be approaching a constant value of ~ 200 MeV 
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2T| , |6l|p6| . The pion mass was extracted from (7r°(x)7r°(0)) and the resulting versus ruf 



showed noticeable curvature for the quark masses used, which were in the range 0.02 — 0.22. 
Therefore, the extrapolation to mj = was done using only the three lightest quark masses: 
0.02, 0.06 and 0.10. Figure |14| updates the earlier graph in with more data at Lg = 48 



and a new point at Lg = 64. The additional data does show a behavior that is more 
consistent with a monotonic decrease of m^(mj = 0) with increasing Lg that than seen in 
our earlier study However, the dependence on L<, shown in Fig. [T^ still cannot be 



described by a single falling exponential and, for large Lg is falling quite slowly. In this 
section we will probe this issue and others related to the chiral limit, using information from 
our simulations at both (3 = 5.7 and 6.0 for many values oi nif and Lg. 

Figure ^ shows results for versus nif for 8'^ x 32 lattices at /3 = 5.7 with Lg = 48, 
including results for rrif = 0.0. The pion mass is extracted from three different correlators: 
(7r"(x)7r"(0)), (v4S(x)Ag(0)) and (7r'^(x)7r'^(0)) + (a(x)a(0))c. For 0.02 < ruf < 0.15, the 
pion masses extracted from the different correlators are in good agreement. As m/ — >■ the 
masses begin to disagree, presumably due to the differing contributions of the topological 
near-zero modes to each correlator. Table pCXl| gives our fitted pion masses for mj < 0.01. 



While the different correlators generally have reasonable values for x^/dof the fitted masses 
disagree substantially. For large enough separation of the interpolating operators, the three 
correlators should give the same mass. However, we cannot take this large separation limit 



in our finite volume. The results in Table pCXl| are the apparent masses as determined from 



fitting to the correlators for finite separation of the interpolating operators. 

The lines drawn in Figure ^ are from correlated linear fits to using rrif = 0.02 to 0.1. 
The dotted line is for ml from (7r'^(x)7r'^(0)), the solid line for {Aq{x)Aq{0)) and the dashed 
line for (7r''(a;)7r'^(0)) + (a(x)(T(0))c. The fit results are 

= 0.053(8) +4.76(12)m/ xV^of = 0.7 ± 2.5 (69) 



m 



0.042(8) +4.90(9)mj x /dof = 0.01 ± 0.25 (70) 
0.037(8) + 5.04(6)m/ xV^of = 1.7 ± 2.8 (71) 



for (7r"(x)7r"(0)), {A^{x)A^{0)) and (7r"(x)7r"(0)) + {a{x)a{0))c respectively. Note for large 
mass, nif > 0.1, (7r"(x)7r"(0)) + ((T(x)cr(0))c gives a mass that is systematically higher than 
that implied by the other two correlators, likely due to contamination from heavy states 
present in (cr(x)(j(0))c. 

Figure |T3| shows the pion effective mass from the three different correlators for 8^ x 
32 lattices at /3 = 5.7 with Lg = 48 and nif = 0.0. Reasonable plateaus are present 
in {Aq[x)Aq{0)) and (7r'^(a;)7r"(0)), although the value for the effective mass is markedly 
different. The (7r"(x)7r"(0)) + {a{x)a{0))c effective mass becomes very small for intermediate 
values of t. Figure ^ is a similar plot, except for rrif = 0.04. Here the effective mass plots 
show nice plateaus and consistent results. This supports the presence of topological near- 
zero modes affecting the various correlators in different ways and provides an example where 
nice plateaus do not assure a correct asymptotic result. 

As a final step in demonstrating the zero mode effects in the various correlators, in 
Figure |13 the evolution of (7r''(x)7r''(0)), {a{x)a{0))c and (Ag(x)y4g(0)) is shown for 8^ x 32 



lattices at P = 5.7 with Lg = 48 and nif = 0.0. These correlators are from a point source 
to a point sink and the zero spatial momentum component is taken for the sink. The sink 
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is at a separation t = 8 from the source. The correlators (7r"(x)7r°(0)) and — (cr(x)(j(0))c 
show very large fluctuations, which are common to both correlators, showing the presence of 
topological near-zero modes. These large fluctuations are clearly dominating the ensemble 
average for the correlators at this separation, t = 8. The {Aq{x)Aq{0)) correlator does not 
show large fluctuations where (7r'^(a;)7r"(0)) and (o"(x)cr(0))c do, making the topological near 
zero mode effects smaller for this correlator, as expected from the theoretical discussion of 
the previous subsection. Figure ITBI is a similar plot, for the same configurations, except with 



rrif = 0.04. There is no evidence for a large role being played by the topological near- zero 
modes. 

Similar results have been obtained for simulations on 16^x32 lattices at /5 = 6.0 with Lg = 
16. These lattices have essentially the same spatial volume, in physical units, as the previous 
8^ X 32, /? = 5.7 lattices since the lattice spacing is half that for P = 5.7. Figure ^ shows ml 
for (Ag(x)Ag(O)), (7r"(a;)7r''(0)) and (7r''(x)7r'^(0)) + (cr(x)a(0))c. For the smallest mj points, 
0.0 and 0.001, all three correlators give different results. For larger values of ruf, the pion 
mass from {Aq{x)Aq(0)) and (7r'*(x)7r°(0)) agree, while the (7r'*(x)7r°(0)) + {a(x)a{0))c pion 
mass is systematically high, likely due to the contribution of the heavy states in {cr{x)a{0))c. 
The lines drawn in Figure |19| are from correlated linear fits to using rrif = 0.01 to 0.04. 
The dotted line is for from (7r'^(x)7r'^(0)), the solid line for {Aq{x)Aq{0)) and the dashed 
line for (7r''(a;)7r'*(0)) + {a{x)a{0))^. The fit resuhs are 

ml = 0.0132(20) + 3.07(6)m/ xV^of = 5.0 ± 5.0 (72) 
ml = 0.0098(20) + 3.14(9)m/ xV^of = 0.03 ± 0.30 (73) 
= 0.0020(26) + 3.56(8)m/ xV^of = 0.06 ± 0.51 (74) 

for (7r"(x)7r''(0)), (Ag(x)Ag(O)) and (7r"(x)7r''(0)) + {a{x)cT{0))^ respectively. 

Figure |20| shows effective mass plots for the pion from the three correlators for mj = 0.001 
and Figure ^ is for mf = 0.01. Both figures show reasonable plateaus, even though there are 
differences in the final fitted masses. We have also studied the evolution of the correlators 
at a fixed t for these /3 = 6.0 lattices and see clear topological near-zero mode effects as were 
seen at /3 = 5.7. 

Thus, investigating the chiral limit of domain wall fermions in quenched QCD by mea- 
suring the pion mass is made difficult by the presence of topological near-zero modes. One 
component in the somewhat large values of ml plotted in Fig. |l^ is the effect of topologi- 
cal near zero modes. As can be seen by a comparison with Table |XV1| the results we find 



from the correlator {Aq{x)Aq{0)) for m^(mj = 0) are about 1 1/2 standard deviations lower 
for Lg = 32 and 48. This is likely a systematic bias caused by the greater influence of the 
topological near-zero modes on the (7r'*(a;)7r"(0)) correlator. Unfortunately, these effects may 
also enter in the other correlators that can give m^^, at least for the source-sink time separa- 
tions currently accessible. With this large distortion due to the topological near-zero modes, 
there we cannot determine the chiral limit by extrapolating to the point where ml van- 
ishes. Subtler finite volume effects and possible quenched chiral logarithms are completely 
overshadowed by the singular nature of the basic quark propagators for small mj. 

In many ways, the presence of these topological near-zero modes is a welcome change 
from other lattice fermion formulations because they are a vital part of the spectrum of any 
continuum Dirac operator. However, in order to further investigate the chiral limit, they 
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must be removed, or at least suppressed. Without adding the fermionic determinant to the 
path integral, we can suppress the effect of topological near-zero modes by going to large 
volumes. 



C. The pion mass for larger volume 



Having seen clear evidence for topological near-zero modes in the measurements of the 
pion mass for lattices with a physical size of ~ 2 Fermi, we have worked on a larger physical 
volume, ~ 4 Fermi, to suppress the effects of these modes. As we saw in Section |^ from 
studying (gg), the effects of the topological near-zero modes were dramatically reduced for 
larger volumes. Here we present results for the pion mass from simulating with 16^ x 32 
lattices at /? = 5.7 and = 48. 

Figure |^ shows plotted against rrif for these runs. In contrast to the smaller volume 
8^ X 32 result shown in Figure ^, all three correlators now give the same results for the pion 
mass, within statistics (Table PCX1| ). The larger volume has clearly reduced the effects of 
the zero modes. Further evidence of the consistency of the mass from the three correlators 



is shown in Figure |22|. Here, for each nif, the average value of m,r is calculated and then 
the deviation from that average, for each correlator, is plotted. For each rrif, {Aq{x)Aq[0)) 
is offset to the left and the (7r"(a;)7r'^(0)) + (cr(x)(T(0))c to the right for clarity. 

Figure |2^ shows the effective mass from each of the three correlators for m/ = 0. In 
contrast to the smaller volume case, the effective masses have quite similar values and lead 
to the same fitted mass, within errors. As a last comparison with the small volume. Figure 
24| shows the three correlators at a time separation of t = 8 as a function of configuration 
number. Little if any effect of topological near-zero modes is seen. Thus, we conclude that 
this larger volume has suppressed these effects as expected. 

Having established that a consistent pion mass can be determined from our fitting range, 
we discuss the result of linear fits of as a function of rnj. We have done correlated linear 
fits of to rrif for each of the correlators, using a variety of different ranges for rrif in the 
fit. The resulting P^r degree of freedom is shown in Figure |25|, including the jackknife 
error on the x^- (The plotted error bars are the ilex errors from the jackknife procedure and 
do not mean that can become negative.) The pion propagator for rrif = 0.0 and 0.04 was 
measured on the same set of configurations, with some of the ruf = 0.08 propagators also 
measured on those configurations. The rrif = 0.02, 0.06 and 0.10 points were all measured 
on the same configurations, along with the remaining 0.08 propagators. Thus, these points 
are less correlated in m/ than the corresponding measurements on the smaller volumes. 

Now let us discuss the quality of these fits. Given the significant upward curvature of 

for rrif > 0.1, seen for example in Figure 0, we limit the mass range to rrif < 0.08. If 
we do not include the lightest masses and fit the points with 0.02 < rrif < 0.08, as shown in 
Figure ^ we obtain acceptable values for per degree of freedom for all three correlators. 
Specifically using the mass range rrif = 0.02 to 0.08, the fits to from the correlators 
(7r'^(x)7r'^(0)), (Ag(x)Ag(O)) and (7r''(x)7r'^(0)) + (a(x)a(0))c are 

ml = 0.044(5) + 4.75(5)m/ xV^of = 1.4 ± 3.6 

ml = 0.051(3) + 4.68(4)m/ xV^of = 1.4 ± 1.4 (75) 
ml = 0.049(3) + 4.70(5)m/ xVdof = 2.5 ± 2.4. 
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However, given our confidence that this larger 16^ volume permits the reliable calculation 
of the pion mass for smaller values of m.^ we can also attempt a linear fit in the entire range 
0.0 < rrif < 0.08. For this mass range, we find 

ml = 0.042(3) + 4.77(3)m/ xV^of = 1.3 ± 3.5 

ml = 0.044(3) + 4.75(5)m/ xV^of = 4.3 ± 2.6 (76) 
ml = 0.042(4) + 4.82(6)m/ xV^of = 4.4 ± 3.0 

for the correlators (7r"(x)7r'*(0)), {A^{x)Aq{0)) and (7r"(a;)7r"(0)) + {a{x)a{0))c respectively. 
The {Aq{x)A^{0)) and (7r"(x)7r"(0)) + {a{x)a{0))c fits suggest that ml is not linear in this 
mass range. While the (7r"(x)7r"(0)) fit is acceptable, as can be seen from a careful examina- 
tion of Figure |^, this acceptable fit comes because the mf = 0.02 point lies somewhat below 
while the rnj = 0.0 lies somewhat above the masses obtained from the other two correlators. 
Since the smaller volume studies suggest that the (7r'^(a;)7r'^(0)) correlator is most sensitive 
to zero modes and such an upturn for small mass is the effect of zero modes seen at smaller 
volume, this could easily be a remaining zero mode distortion. 

It is difficult to draw a firm conclusion from the relatively large correlated x^/dof pre- 
sented in Eq. |76|. As is indicated by the errors shown, these x^/dof are not reliably known. 
However, the comparison of the x^/dof between Eqs. |75| and may be more meaningful. 
We attribute significant weight to the fact that the lightest mf = point lies below the value 
predicted by a linear extrapolation from larger masses as can be easily seen in Figure |^. 

We conclude that a linear fit does not well represent our data over the full mass range 
m/ = 0.0to0.1. Of course, non-linearities for larger masses can come from a variety of 
sources including terms from the naive analytic expansion in powers of mj. However, for 
small mf, linearity is expected for large volumes in full QCD. In contrast, in the quenched 
approximation the absence of the fermion determinant may result in complex and more 
singular infrared behavior. For example, it has been argued that a quenched chiral logarithm 
can appear in ml versus mf for quenched QCD p8|-[rO[1. The results just presented may be 



evidence for some non-linear behavior of this sort. 

Because of the poor linear fits found for small mf, our data does not allow a determination 
of the location of the chiral limit for quenched domain wall fermions by a simple extrapolation 
of ml. Even with the suppression of topological near-zero mode effects that has been 
achieved by going to larger volume, further theoretical input may needed if we are to deduce 
mres from these measurements of ml. In the next section we will discuss our determination 
of the location of the chiral limit using other techniques and then return to the question of 
the behavior of ml with mf. 

VI. THE RESIDUAL MASS 

A. Determining the residual mass 

In this section, we discuss our determination of mjes using the low-momentum identity 
in Eq. |ll[ This can done by calculating the ratio 
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as a function of t (no sum on a), where 7r'^(0) is a source evaluated at t = but possibly 



extended in spatial position. This ratio was first used to determine m^es in Ref. []10| and later 
in Refs. [^,^]. Our results are consistent with this earlier work, but a much more detailed 
study is undertaken here. For t outside some short-distance region, t > tmm, R(t) should be 
simply equal to rrires- Using R{t) for very large t gives rrires as the coupling of the pion to 
the mid-point pseudoscalar density divided by its coupling to the wall pseudoscalar density. 
Of course, rrires is an additive contribution to the effective quark mass at low energies which 
effects all low-energy physics, not just the pion. To understand how large t must be. Figure 
PB| shows a typical good plateau and a poor one. Results are shown for 8^ x 32 lattices 
with rrif = 0.04 and (3 = 5.7 for Lg = 32 and 48. The good plateau is obtained from 335 
configurations for Ls = 48, while the poor plateau is obtained from 184 configurations for 
Ls = 32. The fewer measurements for Ls = 32 likely is the cause for the upturn in the data 
at large t and adding more configurations at this Lg should improve the signal. 

From observing the onset of the plateaus in our data, we calculate mres from the ratio in 
Eq. |73 using the range 4 < t < 16 for /5 = 5.7, 6 < t < 26 for /5 = 5.85, and 2 < t < 16 for 
f3 = 6.0. The jackknife method is used to measure the statistical uncertainty and our mres 
results a.t (3 = 5.7, 5.85 and 6.0 are listed in Tables pCXII| and |XX1I1| . For most data sets, 
nice plateaus can be seen over the selected range, while for the few others with the poor 
plateaus, using a different range could change the results by < 5%. We have also measured 
mres for different values of mj for j3 = 5.7 on 16^ x 32 lattices with Lg = 48. Table PCXIII] 
gives the results and shows that the residual mass has little dependence on the input quark 
mass, reflecting the expected universal character of mres- Our /3 = 6.0 results for mres appear 
to be a consistent extension of the values plotted in Figure 5 of Ref. for = 4, 6 and 
10. 

The Lg dependence of mres is of vital importance to numerical simulations with domain 
wall fermions. Without the effects of topological near-zero modes, quenched chiral logs and 
finite volume, m^(m/ = 0) should be proportional to mres and should vanish with mres as 
Lg —>■ oo. However, in Section we discussed how topological near- zero mode effects alter 
(7r"(a;)7r'^(0)) and can distort the value of m^(m/ = 0) for large Lg shown in Figure |14. 
By measuring the ratio in Eq. |7^, we can determine mres for non-zero rrif and suppress all 
these effects which make the ^ limit problematic. This allows us to study the Lg 
dependence of mres, to which we now turn. 

From the two values of Lg shown in Figure we see that the residual mass for 8^ x 32 
lattices &t f3 = 5.7 falls from 0.0105(2) to 0.00688(13) as Lg is increased from 32 to 48. This 
is in sharp contrast to the almost identical results for m^(mj = 0) at these two values for 
Lg (Figure ITl). The overlap of the surface states is significantly suppressed, as expected, 
even at this relatively strong coupling. We have not pursued the asymptotic behavior for 
large Lg at P = 5.7, due to the large values for Lg required, but instead have studied this 
question for /3 = 6.0. 

Figure ^ shows a similar study of the Lg dependence of the residual mass for 16^ x 32 
lattices with rrif = 0.02 and /3 = 6.0. The number of configurations used is modest for 
the larger values of Lg. We have used the factor Z5(MS,2GeV) = 0.619(25) obtained by a 
combination of non-perturbative renormalization and standard perturbation theory ||72|J73[] 
to convert the plotted values of mres into MeV. The value of mres is decreasing with Lg for 
all values of Lg, but is poorly fit by a simple exponential. In particular, an exponential fit 
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using all values of Lg gives 

TUres = 0.0068(4) exp(-0.094(4)L,) , xVdof(3) = 32 (78) 

which clearly does not match the measured values. Adding a constant to the fit gives 

m,es = 0.00032(3) + 0.018(3) exp(-0.181(13)L,), xVdof(2) = 4.1 (79) 

where again all values for Ls were used. Even if this is the correct asymptotic form, the 
value of rrij-es for — oo is very small, 1 MeV. 

We have also tried fitting the largest three Ls points to a simple exponential and find 

m,es = 0.0012(2) exp(-0.032(6)L,), xVdof(l) = 0.074. (80) 

Our data is consistent with the residual mixing vanishing exponentially as — > oo, but the 
0.032 coefficient in the exponent of Eq. ^ is quite small. Of course, we can easily obtain an 
excellent fit to our five points if we include a second exponential. For example, as shown in 
the figure, the five points fit well to two-exponential function 

TUres = 0.038(16) exp(-0.26(4)L,) + 0.0010(3) exp(-0.027(7)L,), xVdof(l) = 0.1. (81) 

Our measurements do not demonstrate a precise asymptotic form for function 
of Lg. However, we do see mres decreasing for large Ls until, for Ls > 24, it has become 
so small as to be essentially negligible for current numerical work. For jS = 5.7 at Ls = 
48, mres is 0.074(5) in units of the strange quark mass, while for (3 = 6.0 at Ls = 16 it is 
0.033(3)ms. In the latter case, where we know the renormalization factors, mres in the MS 
scheme at 2 GeV is 3.87(16) MeV. Thus, even though more simulations will be needed to 
get the precise asymptotic form, we find domain wall fermions having the expected chiral 
properties for large Ls, even for lattice spacings of around 1 Fermi. 

In the next subsection, we will use the values of mres that have just been determined to 
investigate further the rrif dependence of m^, looking in particular at possible non-linear 
behavior as m^ + mres 0. Here we would like to discuss a simpler consistency check on the 
values of mres just obtained. For the tt, p and nucleon we have established good linear mj 
behavior for larger values of mj with slopes and intercepts given in Tables PCV1| and [XVI 1| . If 



the only effect on these masses of changing Ls is to change the effective quark mass through 
the corresponding change in mres, then we should be able to relate the the differences in the 
intercepts given in these tables to the product of the corresponding slope times the change 
in mres given in Tables pCXll| and pCXlli| . 



While this comparison shows no inconsistencies, the errors in the intercepts are typically 
too large to permit a detailed confirmation. For example, the difference in intercepts for 
m^ at /? = 6.0 between Lg = 16 and 24 is 0.0004(30) while the difference predicted from 
the slope and the measured change in mres is 0.0020(2). The best test of this sort can be 
made using the actual value for m^r determined at /5 = 5.7 and m^ = 0.04 for Lg = 32 and 
48. Here the difference of the masses squared is 0.012(6) while the prediction from the slope 
and change in mres is 0.0176(12). Thus, we can demonstrate consistency with the expected 
behavior but cannot make a definitive test. 
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B. The residual mass and versus mj 



The definition of mj-es and its measurement mean tliat we fiave determined tlie value of 
rrif for which the pion should become massless if the domain wall method is successfully 
representing the chiral limit of the underlying theory. We can now return to the question of 
the dependence of on ruf, starting with the 16^ x 32 simulations at /? = 5.7 and Lg = 48. 
Recalling Figure ^, we found that the larger volume gave consistent pion mass measurements 
from the three correlators, but was not well fit as a linear function of for two of the 
correlators if the frif = 0.0 point was included. In Figure ^, we have included the value 
of m^es (the starred point) as measured from Eq. |7^. (Its error bar on the horizontal axis 
is a vertical line on this scale.) The solid line is the fit to the {Aq{x)Aq{0)) correlator for 
rrif = 0.02 to 0.08 given in Eq. ^ while the dotted line is for the (7r"(a;)7r'^(0)) correlator for 
rrif = 0.0 to 0.08 as given in Eq. |7^. Thus, we see that linear fits poorly represent the data 
when the mj = 0.0 point is included for the {Aq{x)Aq{0)) and (7r"(a;)7r"(0)) + {a{x)a{0))c 
case and fail for all three correlators when the pion mass is required to vanish at rrif = —m^^s- 

We can make this conclusion more quantitative by comparing our accurate value for 
rrires = 0.00688(13) at = 48 determined on an 8^ X 32 lattice with the naive linear 
extrapolation of m^{mf) to the point = 0. Using the most reliable linear fits obtained by 
excluding the rrif = point in Eq. we obtain the x-intercept values shown in Table |XXIV| : 
-0.0092(12), -0.0108(7) and -0.0104(7) for the (7r"(2;)7r'^(0)), {A^{x)A^{0)) and (7r''(x)7r'^(0)) + 
((j(x)cr(0))c correlators respectively. These differ from this value of m^es by ~ 50% and 2, 5 
and 6 standard deviations respectively. We conclude that there is a significant discrepancy 
between the m/-dependence of these m^r results and the hypothesis that m^(— mres) = 0. 
However, notice that if the rrif = points are included in the linear fits, and the less accurate 
rrij-es from the same volume is used, this discrepancy can be reduced. For example, a linear 
fit to the data from the (7r"(2;)7r'*(0)) correlator in Eq. ^ has an intercept at -0.0088(5) while 
—rUres = —0.0072(9) on the same volume. We believe that such an interpretation should be 
discounted as failing to exploit all the available information. 

Is this significant discrepancy caused by essential non-linearities in the quenched approx- 
imation or by a breakdown of the domain wall method, for example, large O(a^) effects? 
We can address this question by making a similar comparison for (3 = 6.0 where O(a^) 
effects should be significantly reduced. Since we have not investigated a large volume at 
this weaker coupling, we propose to examine the {Aq{x)Aq{0)) correlator because reduced 
zero- mode effects were seen for this correlator in our j3 = 5.7 studies. Using the three light- 
est masses we find x-intercepts of -0.0031(7) and -0.0030(9) for the Lg = 16 and 24 cases 
respectively. Again, these are dramatically farther from the origin than the corresponding 
values of m^^s = 0.00124(5) and 0.00059(4). These are each three standard deviation effects. 
However, they are obtained on independent configurations and together can be viewed as a 
6 standard deviation discrepancy. Thus, if possible finite-volume difficulties are ignored, we 
have again strong evidence for a discrepancy. Rather than decreasing by a factor of four as 
would be expected from an O(a^) error, this fractional discrepancy is substantially larger 
in this /5 = 6.0 comparison. Thus, it is natural to conclude that domain wall fermions are 
accurately representing the chiral behavior of quenched QCD. 

At the beginning of Section |V| we listed possible systematic effects influencing the chiral 
limit for m^. With a measurement of mres we have quantified the role of finite and with 



31 



the larger volume used for /3 = 5.7 we have reduced, if not eliminated the topological near- 
zero modes. We should also have minimized other finite volume distortions of the density of 
eigenvalues, which also influence the pion mass. Finally with the comparison above, we have 
examined the possibility of O(a^) errors. Thus, we now address the question of quenched 
chiral logarithms. Predictions of this particular pathology of quenched simulations were 
made some time ago. There is certainly much data indicating possible support for the 
predictions, but there is disagreement about its conclusiveness, see for example refs. ||7^-|75|. 



Since many other effects must be removed before these subtle logarithms are convincingly 
seen, it is a challenging numerical issue. 

The natural first place to look for quenched chiral logarithm effects is in m^r, but this is 
difficult for Wilson fermions, where the chiral point is not crisply defined for finite lattice 
spacing. For staggered fermions, where the chiral limit occurs when the input quark mass 
is zero, the issue is complicated by the presence of only a single Goldstone pion. In some 
respects, domain wall fermions are an ideal place to look for these effects, except that the 
statistical resolution needed is difficult to achieve with the additional computational load of 
the fifth dimension. In addition, the topological near- zero modes are a much larger quenched 
pathology at moderate volumes. 

As one way of probing the non-linearity in versus mj, we have fitted our data for 
for 16^ X 32 lattices at (3 = 5.7 and = 48 to the form 



= ao(m/ + ai)(l + a2 ln(mj + ai)) (82) 
and the results are given in Table pCXIV|. The fit yields a value for the residual mass (the 



parameter ai above) and the results are quite close to those measured from the ratio of 
Eq. ^ Figure ^ shows the result from fitting (Ag(x)Ag(O)) for mj = 0.0 to 0.08 to the 
quenched chiral logarithm form given in Eq. |8^. We have excluded the larger values of m/ 
from our fits, since higher order terms are needed in Eq. to accommodate the upward 
curvature of our data. While the x^/dof for the logarithmic fit is only marginally better 
than those obtained for the simple linear fits described earlier in Eq. ^ for this same mass 
range, the ability of the logarithmic fit to predict the appropriate rrires value is significant. 

For the simulations at smaller physical volumes, 8^ x 32 at /5 = 5.7 and 16^ x 32 at /3 = 6.0, 
the values for mres measured from Eq. are generally smaller than the x intercepts for the 
linear fits shown in Figures ^ and |19|. This indicates curvature in the direction given by a 
chiral logarithm, but the other phenomena that may be affecting these chiral limits make 
quantitative analysis ambiguous. We note that from (7r"(a;)7r'^(0)) + {a{x)a{0))c seems 
to smoothly curve towards the value of mres from the previous subsection. However, we are 
not sufficiently certain of the absence of zero mode effects in the (7r"(x)7r°(0)) + {a{x)a{0))c 
correlator to describe a logarithmic fit to these cases. 

This nice agreement between the values of mres determined from the location of the 
m^ = point in these fits and that computed by other means earlier in the paper implies 
consistency between our results and the logarithmic form of Eq. Of course, other non- 
linear terms could be used to explain this curvature and, given our statistics, would provide 
an equally consistent description of our data. 

However, our most important conclusion is not related to quenched chiral logarithms, 
but rather to having seen all the expected properties for the chiral limit with domain wall 
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fermions. Once the topological near-zero mode effects are reduced or eliminated, consistent 
pion masses can be measured. A precise measurement of mres is consistent with our 
versus ruf dependence if, for example, a chiral logarithm term is included. In short, domain 
wall fermions are showing sensible chiral properties, even on lattices with a lattice spacing 
of ~ 1 Fermi. 

We have chosen not to pursue an additional method of determining rrires that has been 
proposed in two of our previous publications [^|T^. In that method, one examines the 



integrated Ward-Takahashi identity in Eq. ^ and uses the location of the pion pole in Xmr to 
determine mres- While this technique should be reliable for dynamical fermion calculations, 
e.g. as used in Ref [Q, it does not explicitly allow for the effects of topological near- zero 
modes or possible non- linear behavior of m^(m/) that we have found to be important in the 
quenched approximation. Thus, even though this method gave a result for /3 = 5.7 quite 
close ||2^ to the Lg = 48 value mres = 0.00688(13) presented in this paper, more analysis is 



needed to adequately justify its use in this quenched case. 



C. Eigenvalue properties and rrires 

A comparison of the approximate form of the Banks-Casher relation for domain wall 
fermions given in Eq. ^ with the usual 4-dimensional expression in Eq. suggests a close 
relationship between the parameter 6mi deduced from the i*'* eigenvalue A^^ j of Dh and the 
residual mass m^es- In this section we will explore this relation further making use of an 
exploratory study of the low-lying spectrum of • 



These eigenvalues were calculated for 32 configurations obtained at (3 = 6.0 on a 16^ 
lattice with Lg = 16 listed in Table |I| and beginning with an equilibrated configuration 



from an earlier run. We used the Kalkreuter-Simma I^Sj] method to find the 19 lowest 
eigenvalues on each configuration^. We apply this method to the positive matrix D'jj, and 
then determine the eigenvectors and eigenvalues of Dh by a final explicit diagonalization of 
Dh in the subspace of the eigenvectors of Djj just determined. The details of our application 
of this method and a more complete description of these results will be presented in a later 
publication |]80[ . 



While this method determines both the eigenvalues and eigenvectors, we have chosen to 
examine only the s-dependent, four-dimensional inner products: 

^R/L{s)i,j = E ^^(a^' s)a„,^Pr/l'^{x, s)a^,,, (83) 

X 

where the indices i,j run over all of the 19 eigenvalues while Pr and Pl are the left and 
right spin projection operators defined above Eq. ^ In order to be able to make use of 
the mass dependence of the eigenvalues, we have repeated the calculation of Aj(m/) and 
^R/L{s,Tnf)ij five times on each configuration for the five different mass values frif = 0.0, 
0.0025, 0.005, 0.0075 and 0.001. 



^We thank Robert Edwards whose program formed the basis of the code used in this part of the 
calculation. 
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Here we will describe some of the overall features of this calculation and then examine 
more closely the relation between the parameters 5mj and the value of m^es determined 
earlier in this paper. First we examine the diagonal elements of the matrix T{s) 

^{s), = rR{s),,i + rL{s)i,i. (84) 

This is the contribution to the norm of the 5-dimensional wave function from the 4- 
dimensional hyperplane with a specific value of s. For these low lying eigenvalues, we 
expect that this norm should be concentrated on the s = and s = Lg — 1 walls, which we 
find to be true to good accuracy. For the entire group of 32 x 19 x 5 = 3, 040 eigenvectors 
computed, the ratio of the sum of the norm on the two walls to the minimum value of this 
norm between the walls was always greater than 34, J\f{0) +M{Ls — 1) > 34 ■ Af{smin)- The 
median value for this ratio was 744. Thus, the general framework upon which the domain 
wall formalism rests appears approximately valid. 

As a test of our method for determining the eigenvectors, we evaluate the left- and 
right-hand side of the symmetry relation, Eq. between pairs of eigenvectors on a given 
configuration. The resulting equality: 

iAH, + AH,,){AHA^5\AH,,) = (AH,|(2m;Q(-) +2Q(-P))|A^,,) (85) 



provides a good test of our diagonalization procedure. The vectors needed to evaluate 

this expression are eigenvectors of the Dirac operator Dh, not Djj. We determine the 
eigenvectors of Dh by diagonalization within the 19-dimensional subspace found by applying 
the Kalkreuter-Simma method to Djj. In the event that the 19*'* and 20*'* eigenvalues of 
Djj are nearly degenerate (not entirely unlikely given the expectation that the eigenvalues 
of Dh occur in ±Ah pairs), this truncated, 19-dimensional subspace will not be spanned by 
eigenvectors of Dh- It will contain 18 valid eigenvectors and a 19*'* vector, orthogonal to the 
rest but not an eigenvector of Dh- This "spurious" eigenvector can be reliably removed since 
it will give an "eigenvalue" whose square does not agree with any found for Djj. We remove 
such eigenvectors from our test of Eq. 85 and, for uniformity, the 19*'* eigenvector in the case 
that no spurious eigenvector occurs. There are then 32 x 5 = 160 instances where we can 
check 18^ independent elements of Eq. We find that 95% of these 51,840 comparisons 
have a fractional error below 5%. The few cases with significantly worse agreement, result 
from infrequent near degeneracies which challenge the Rayleigh-Ritz method on which the 
Kalkreuter-Simma algorithm is based. 

For most configurations there are easily identified zero modes. Typically the few lowest 
eigenvalues have eigenvectors all of which are bound to the same wall, either s = or 
s = Ls — 1- The corresponding matrix elements {AH,i\^5\AH,i) all have the same sign and 
are within a few percent of 1, showing precisely the structure expected in a four-dimensional 
theory as summarized in Eq. |3^. 

The potential of the domain wall method is nicely displayed by examining the properties 
of one of our better configurations. In Fig we show the magnitude of the elements 
of the matrix {AH^iir^lAHj) in a three-dimensional plot. Note the five zero-modes in this 
configuration are easily recognized. Each has a diagonal matrix element of within 1.5% 
of 1 and matrix elements with other vectors all of magnitude below 0.06. The values of Aj 
for these five eigenvalues all lie in magnitude below 3.6 x 10~^ while the remaining paired 
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eigenvalues lie between 0.028 and 0.093. In Figures ^ and ^ we show the s-dependence of 
the first two zero-modes and the first pair of non-zero eigenvectors, numbers 5 and 6. One 
sees precisely the expected behavior. Both zero-modes are bound to the same wall (as are 
the other three zero modes) while the two paired non-zero modes are nearly symmetrical 
between right and left. This is clearly identified as a configuration with topological charge 
V = +5. 

Of direct interest in this section is the mass dependence of A{mf)H,i and a quadratic fit of 
the sort proposed in Eq. For the small masses we have used, this quadratic form provides 
an excellent fit, after some re-sorting of eigenvalues is performed to account for infrequent 
level crossings as is varied. In order to avoid the possibility that these level crossings may 
have pushed a needed eigenvalue up to beyond number 19, we have excluded those quadratic 
fits which contain the largest eigenvalue at m/ = for each of the 32 configurations. The 
resulting root-mean-square of the fractional differences between the left- and right-hand 
sides of Eq. P is very small. The average root-mean-square of the fractional difference is 
1.3 X 10"^ while the largest value is 4.3 x 10~^. 

In Figure |31| we present a histogram of the distribution of fit parameters for the 18 x 32 = 
576, Srrii values that we obtain. The majority of values are quite small, very much on the 
order of mres- While a few larger values of 6mi are seen (the largest is 0.0660), the median 
of the distribution is Sm = 0.00147 which is remarkably close to the value of m^es = 0.00124 
found earlier for this value of (3 and Lg. 

The 4-dimensional expression for (ipip) in Eq. |3^ and the 5-dimensional result in Eq. ^ 
as a function of m = Zruf must agree in the continuum limit after a rescaling and overall 
subtraction. This must be true even if mj-es is held fixed in physical units as a — 0. Therefore, 
in the limit of zero lattice spacing, the histogram shown in Figure 31 must approach a delta 
function so that 6mi has the unique value mj-es- Thus, we might interpret the width of 
the distribution in Figure ^ as a result of O(a^) effects. The large size of the fluctuations 
relative to the central value is presumably a result of the small central value produced by 
our quite large separation of 16 between the walls. 



VII. HADRONIC OBSERVABLES 



We can now use the results of the previous sections to compute a variety of hadronic 
properties. In this section we will discuss two topics: the evaluation of the pion decay 
constant f^, and the scaling properties of the nucleon to p mass ratio. The first topic is 
of greatest interest since we can compute the pion decay constant using two independent 
methods, one of which depends directly on the residual mass determined in Section |V1|. The 
close agreement between these two approaches provides a very important consistency check 
of the analysis and results presented in this paper. 



A. Calculation of /, 



In the conventional continuum formulation, the pion decay constant /tt is defined through 
the equation 
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(0^^7W>^'(P)) = U r^"'" (86) 
^ ^ '2EJp) 



where the fields ip and ip are interpreted as conventional, Hilbert space quark operators and 
the pion state obeys the non-covariant normalization {7t{p')\7t{p)) = 5^{'p — p'). To be con- 
crete we adopt the Minkowski metric g^'^ with signature (—1, +1, +1, +1) and a Minkowski 
gamma matrix convention in which 7°^ is anti-hermitian and {7m;7m} = 2(7^*^. With this 
normalization, f^^ k, 130 MeV. 

Following the usual methods of lattice gauge theory, we evaluate matrix elements of the 
two-quark operator appearing in Eq. |86| with the Euclidean time dependence resulting from 
use of the evolution operator e"^* where H is the QCD Hamiltonian. Thus, we choose to 
evaluate 

Tr{e-^(^-*) /c?3xzV^70,75rV'(f,t)e-^*2V570,75rV'(0,0)} 
^'^^ = T^T. ^ TWiF^} 

where no sum over the fiavor index a is intended and the time t is assumed sufficiently large 
that only the pion intermediate state contributes. 

The continuum operators in Eq. ^ are easily represented as lattice, Euclidean-space ex- 
pressions once the usual transition to a Euclidean-space path integral has been performed. In 
particular, the operators ip{x) and %l){x) = ijj^{x)''y'^ are replaced by the Grassmann variables 
q{x,t) and q{x,t) respectively. Thus, we extract from the usual Euchdean correlation 
function: 



-e""'^* = {Jd^x g7°7^rg(f , t) q^^^H^qiQ, 0)) (88) 



where now Euclidean gamma matrices appear, obeying {7'^,7'^} = 25^^. Here we have 
introduced the Grassmann variables q and g defined earlier in this paper so the axial current 
appearing in Eq. ^ is explicitly constructed from the five-dimensional quark fields \l/ and 
\1/ restricted to the s = and s = — 1 walls. This "local" current, is not conserved in 
the full five-dimensional theory so the factor Za appearing on the left hand side of Eq. ^ 
is needed to make a connection to the continuum axial current. 

The conserved current defined in Eq. ^ must approach the corresponding, partially 
conserved continuum current with unit normalization, when the continuum limit is taken. 
Thus, to order a^, the low energy matrix elements of and must be proportional: 
A""^ = ZaAIj^. While we have computed using the local current we have also compared 
that current to the partially conserved domain wall axial current allowing an accurate 
determination Z^- 

In addition to the procedure just described, there is a second, independent method that 
we have used to compute /^r. Here we use the Ward-Takahashi identity to relate A^^ and the 
pseudo-scalar density J5 : 

A^A'^^ix) ^ 2(m; + m,,,)J^{x) (89) 

an expression valid for low energy matrix amplitudes . In particular, we have replaced the 
usual midpoint term in the exact identity of Eq. |^ by its low energy limit: 2mTesJ^{x). Thus, 
we can also obtain from the correlation function: 
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J 7r 



3 



- ^ Jd^x q^H^'qix, t) qjH^'qiO, 0)) (90) 



where again no sum over the flavor index a is intended. This formula involves no renormal- 
ization factors but requires knowledge of the residual mass rrij-es induced by mixing between 
the walls. Thus, a comparison of the values for obtained from Eqs. |8^ and ^ provides a 
critical test of the analysis presented in this paper. 

We will now discuss these two calculations of in detail. To measure the value for the 
renormalization factor Za, we compare the amplitudes of two-point functions C{t) and L{t) 
defined as 

C{t + 1/2) = J2{A^,{x,t) 7r"(0,0)) 

L{t) = Y.{A^oi^,t)7r^{6,0)). (91) 

X 

The 1/2 in the argument of C(t + 1/2) in Eq. |^ comes from the fact the conserved axial 
current Al^{x) is not the current at lattice site x but instead the current carried by the link 
between x and x + fi. We take appropriate arithmetic averages to solve the problem that 
C(t + 1/2) and L{t) are not at the same location. To avoid as much systematic error as 
possible, we define Z^it) as 



C{t + 1/2) + C{t - 1/2) , 2 C{t + 1/2 
-"^ ^ 2\ 2 Lit) ^L{t) + L{t' 




(92) 



For t ^ a^^, C{t)/L{t) behaves like a constant which can be identified with Za- Both terms 
in Eq. ^ estimate this value without 0{a) error. The average of these two, incorporated in 
Eq. further eliminates a portion of the O(a^) error. 

Figure ^ shows the ratio Z^(t) defined in Eq. ^for both a 16^ x 32 lattice with Lg = 16, 
and (3 = 6.0 as well as the same quantity for a 8^ x 32 lattice with Lg = 48 and (3 = 5.7. We 
determine the value for the renormalization factor Za by calculating the average over two 
ranges of t: 4 < t < 14 and 18 < t < 28, chosen to avoid the largest time separation t ~ 16 
where the errors are quite large. A jackknife error is determined, to compensate for possible 
correlation between the numerator and denominator in Eq. p2| . 

The results for Za ai (3 = 6.0, 16'^ x 32, M5 = 1.8 and with different values of Lg are 
listed in Table pCXV|. The data shows little Lg dependence, as should be expected. Figure 



32| , also shows our result of Za = 0.7732(14) found for the 8^ x 32 lattice with /3 = 5.7, 
Lg = 48, M5 = 1.65, nif = 0.02. 

The results for the amplitudes for the axial vector current correlator and the pseudoscalar 



density correlator at /? = 5.7 and 6.0 are given in Tables pCXVI[ - ^XX| . They are obtained 



from the point-source correlators using a conventional 2-parameter fit with the pion masses 
extracted concurrently. We also list in the same tables the results for f^^ as a function of 
determined from the corresponding correlators with the help of Za and rrires (Tables pCXli| - 



XXIII|) . These values of f^, have been converted to physical units using the measured p mass 
discussed in Section |TT|, extrapolated to the chiral limit mf + m^es = 0. 

Next, we use a linear fit in to evaluate f.^ for two values oi mf. To obtain a value of 
close to that for the physical pion, we go to the chiral limit ruf + nires = 0. For fx we 
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choose for rrif that value which gives m.^/mp = 0.645. In determining /,r for the physical 
pion state, we did not attempt to use a value of m/ giving the physical value for the ratio 
m^/mp = 0.18 since we do not adequately know the ruf dependence of this ratio in the 
relevant region. These linear fit parameters as well as the resulting values for /^^ and fx are 
summarized in Table pCXXl| . The errors given in the tables are obtained from the jackknife 
method. 

Figure ^ shows the values for /^r at /5 = 5.7, 8'^ x 32, Lg = 48 as a function of m/ 
and the linear fits through all the mj points. The results obtained from the pseudoscalar- 
pseudoscalar correlator are higher than those from the axial-axial correlator. The two linear 
fits give U = 127(4) MeV, fx = 145(4) MeV and U = 132(4) MeV, fx = 154(4) MeV 
respectively. When the lattice volume is increased to 16^ x 32 (Figure ^), the difference 
between the linear fits from the two methods becomes smaller. We obtain f.^ = 133(4) MeV, 
fx = 149(2) MeV and U = 125(4) MeV, fx = 149(2) MeV from the two correlators. The 
values for fT^{'mf) obtained from the two methods should agree for all values of since 
they are related by a Ward-Takahashi identity that should become exact in the continuum 
limit. Presumably the visibly different slopes seen in Figures ^ and |3^ are the result of 
order a? errors. 

We also calculate f^, at a weaker coupling. Figure ^ shows our results for /5 = 6.0, 16^ x 32, 
Lg = 16 on 85 configurations. The two independent calculations give very consistent results. 
We have /^r = 137(11) MeV, fx = 156(8) MeV from the axial vector current correlator 
and almost the same values from the pseudoscalar correlator. Our results for f.^ at both 
f3 = 5.7 and 6.0 agree well with the experimental value of ~ 130 MeV, while the values for 
fx may be somewhat smaller than the experimental value of ~ 160 MeV as is expected from 



quenched chiral perturbation theory arguments [B9] and naive scaling considerations [81 



Note, in Table pCXXi] we also list fxl f-w with jackknifed errors for the ratio. Here the 
statistical errors are now well below the systematic errors that might be expected in the 
m/ + mres — > extrapolation. The values shown for fx/ fn agree on the 5% level between 
methods of determination and different lattice spacings but are systematically below the 
experimental value of 1.21. 

This same analysis was done using the amplitudes calculated from the point-source cor- 
relators but making a 1-parameter fit using the pre- determined pion masses computed from 
the more mass accurate measurements based on the wall-source correlators. This method 
gives consistent results with slightly smaller errors. The results are not listed here. 

The reasonable agreement of our domain wall results with the experimental values and 
their relative insensitivity to a is encouraging. Similar results were obtained at /3 = 6.0 for 
smaller values of with somewhat larger errors in Ref. PDI. Of special interest here is the 



comparison that we make between the two methods of determining /^r, which is done here 
for the first time. As can be seen from Eq. |9^, the determination of f-j^ from (7r"(x)7r'*(0)) 
depends directly on mres- Thus, the comparison of these two methods is an important check 
of our understanding of the chiral properties of the domain wall formulation. The ratio of 
these two quantities extrapolated to the point mf + m.^cs = provides an interesting figure of 
merit for the present calculation. We find {fn)?? I {,f-K)AA = 1.00(10) and 0.96(10) for Ls = 16 
and 24 respectively. However, if instead of the values of mres given in Table pCXll| , we use the 
x-intercepts -0.0031(7) and -0.0030(9) quoted earlier and obtained from the {Aq{x)Aq{0)) 
values of m^, we find {f-K)ppl W-k)aa = 1.20(12) for both the = 16 and 24 cases. While 
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these ratios each differ from 1 by two standard deviations, they are independent calculations 
and demonstrate the good chiral properties of domain wall fermions. 



B. Continuum limit of mj^jmp 



Here we combine the hadron mass results tabulated in Section |II| to examine the behavior 
of the nucleon to p mass ratio as (5 varies between 5.7 to 6.0. First we evaluate and nip 
in the limit w^/ + vn^cii = 0. We did not use the value oi nif which gives the physical ratio, 
mTr/rrtp = 0.18 for the reasons outlined in the previous section. In Table pCXXli] we give the 



resulting mass ratios as well as the lattice spacings in physical units as determined from nip 
evaluated at m/ + m^es = 0. Note, no contribution to the quoted error for these mass ratios 
arising from the uncertainty in this choice of has been included. 

The relatively large variation of m^r/mp with (3 suggests that the errors shown in Ta- 
ble pCXXll] may be underestimated and makes a simple o? extrapolation to the continuum 



limit somewhat uncertain. Nevertheless the result of such an extrapolation to a —> is 
nifqjnip = 1.37(5). Perhaps more interesting is a comparison with similar quantities com- 
puted at comparable lattice spacings and volumes using Wilson and staggered fermions. For 
staggered fermions at /3 = 6.0 on comparable volumes, one finds [§3'i31 ^^n/^p = 1-47(3), 
a somewhat larger and less physical value than the 1.42(4) and 1.38(4) results obtained here 
for Lg = 16 and 24. However, this comparison is made somewhat ambiguous by the signifi- 
cant finite size effects seen in staggered calculations when going from our 16^ x 32 to larger 
volumes |5^. For Wilson fermions, as reported in Ref. one deduces niN/nip = 1.37(2) 



by linear interpolation between the P = 5.93 and 6.17 values presented, a number remark- 
ably close to our domain wall value. When comparing these values, it is important to recall 
that our 8^ and 16^ spatial volumes are not yet infinite and, as discussed in Section |T|, 
corrections on the order of a few percent are expected. 



C. Determining the chiral condensate (qq) 

Finally we use the results presented earlier to estimate the size of the chiral condensate 
(qq). Naively, one might expect that a physical value for (qq) could be easily identified in 
Table PCX| as the mj-independent term a^, defined in Eq. This quantity represents a 
simple extrapolation of {qq){nif) from large mass down to the point ruf = 0. Given the 
volume independence seen for the parameter oq when comparing the (3 = 5.7, 8^ and 16^ 
volumes in Table pCX|, it is natural to expect that such a choice minimizes the sensitivity to 



the finite- volume zero mode effects that give rise to the more singular a_i term. 

However, there are other issues that must be addressed. Perhaps most obvious is the fact 
that the point my = is not the physical chiral limit because the effects of m^es have been 
ignored. This is easily remedied by using the slope Oi, to extrapolate to the physical point 
nif + ni,.es = 0. The resulting estimate of (qq), in lattice units, is given as the fourth column 



in Table pCXXIIIj . However, because {qq) is a quadratically divergent quantity, we cannot 
expect that all the chiral symmetry breaking effects of domain wall mixing are removed 
by this choice of m/. In contrast to many physical quantities, (qq) receives contributions 
from energy scales much larger than those for which m^es represents the complete effect of 
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chiral symmetry breaking. Thus, we should expect additional contributions to (gg) of order 
^-olLs ~ m^e^ja}. This is born out in Table XXXIIl where we see that the differences 
between (gg) for the two different values of at a given /3 are of the same order as the 
difference between the values with and without the extrapolation to = — mj-es- 

This unwanted ~ my-^^ja^ contribution to (gg) can only be controlled by explicitly taking 
the limit ^ oo. We do not at present have the numerical results to permit such an 
extrapolation. Therefore, we will use the (5 = 6.0, Lg = 24 result as our best approximation 
to such a limit and interpret the difference between the Lg = IQ and 24 values as an estimate 
of the systematic error, ^ 10%. Given the value of Z5(MS,2GeV) = 0.619(25) for f3 = 6.0 
quoted earlier and the results for the lattice spacing in physical units in Table pCXXll] , we 
can determine (gg) in physical units. The results for Lg = 16 and 24, (245(7)MeV)^ and 
(256(8)MeV)^, are included in Table pCXXlll| , where only the statistical error is displayed. 
The agreement between these numbers and phenomenological estimates of the chiral con- 
densate is satisfactory, for example the value of |(MU + rf(i)j^ ^^^y = (229±9MeV)^ obtained 
in Ref. [^. Note the e~°'^'' /a^ uncertainty present in our calculation does not have an ana- 
logue in the properly regulated continuum theory. While (ip'ip) does contain a quadratically 
divergent piece in the continuum theory, this is eliminated for the chirally symmetric choice 
'^quark = 0. This choicc is not available in a domain wall fermion calculation without taking 
the Lg ^ oo limit. Of course, the other lattice methods for directly computing (ipip) have 
equal or more severe difficulties. 

Finally it is interesting to compare the (3 = 5.7 and (3 = 6.0 results for (gg). Since 
we do not at present have a reliable determination of the needed renormalization con- 
stant, Zs, for the stronger /3 = 5.7 coupling, we do not attempt to quote a physical 
value. However, the ratio of the unrenormalized lattice numbers given in Table PCXXIII] for 
(^Q'){Li,=32,/3=5.7)/(^Q')(Ls=24,/3=6.o) = 4.8(2) is reasonably consistent with the ratio expected 
from naive scaling a?L^=32,/3=5.7)/«?L,=24,/3=^o) = 7-4(4). 

Given the values now determined for (gg), /,r and quark mass, it is natural to test the 
degree to which the Gell-Mann-Oakes-Renner relation 

is obeyed. However, the form of this equation reveals an important difficulty. At what value 
of TTif should the ratio m\/{mf + mres) be computed? In full QCD, this ratio becomes a 
constant for small quark mass. As we have seen earlier, this is not the case in the quenched 
approximation where one expects non-linearities. 

We might try to determine the proper treatment of these non-linearities by returning 
to the underlying equation, Eq. ^ from which the Gell-Mann-Oakes-Renner relation is 
derived. However, this is somewhat complex. Both sides of this original equation have 
a mass dependence which comes from the contribution of the pion pole term and other 
physical states, all influenced by the quenched approximation, as well as the quadratically 
divergent terms in (gg) and the contact term in Xi^-n- Thus, while the underlying Eq. ^ 
will be obeyed exactly in our calculation, there is considerable ambiguity in deciding how 
to extract a quenched generalization of the Gell-Mann-Oakes-Renner relation, Eq. p3| . 

Here we will simply compare the right- and left-hand-sides of Eq. ^ by replacing the 
ratio m^/(m/ + mres) by the slope h obtained at larger masses, mj > 0.01 and given in 
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Tables pCVI| , pCVII| and pCIX| . The results from the left hand side of Eq. ^ are given in 
Table XXXIIl . Given our uncertainty in determining (gg) and the significant non-linearities 
we see in m^, the agreement seen between the fourth and fifth columns in Table pCXXlll| is 
within our errors. 



VIII. CONCLUSIONS 

We have presented the results of detailed studies of quenched lattice QCD using domain 
wall fermions, with particular attention paid to the lowest order chiral symmetry breaking 
effects of finite Ls and the behavior of the theory for small values of rrif. A major difficulty 
in studying the small ruf behavior of the theory is the presence of topological near- zero 
modes which are unsuppressed in the quenched theory. These are a result of the improved 
character of the domain wall fermion operator, which has an Atiyah-Singer index at finite 
lattice spacing and oo. However, these zero-modes comphcate the quenched theory 

and demonstrate that the quenched approximation is considerably more treacherous than 
might have been originally expected. We have seen how these modes produce the expected 
1/mj divergence in (qq) for small rUf and distort correlation functions used to measure the 
properties of the pion. By working on larger volumes, we found that the effects of these 
modes were dramatically reduced, as expected. We were then able to see a common pion 
mass determined from different correlators. 

We have determined or constrained the value for the residual mass, rrires, which enters 
the effective quark mass for low-energy physics as rries = rrif + mres, a number of ways 
and found good agreement. The residual mass was measured from the extra, finite Ls term 
in the divergence of the conserved axial current and from the explicitly determined lowest 
eigenvalues of the hermitian domain wall fermion operator. These two determinations agree 
within errors. We have also determined the difference in mres ^ot two values of Lg from 
the pion mass and find this agrees with the results from our explicitly calculated mres- 
Lastly, agreement for as calculated from axial vector and pseudoscalar correlators requires 
knowledge of mres and the agreement serves as a further check. 

While our data for weaker couplings does not clearly demonstrate that mres — 0, we 
have seen it fall to 1 MeV for = 48 at /5 = 6.0. For Ls = 16, a practical value for 
studies of low energy hadronic physics and matrix elements, mres has a value of 3.87(16) 
MeV, roughly 1/30 of the strange quark mass. Even at stronger couplings, where the lattice 
spacing is ~ 1 GeV, we have measured mres to also be about 1/14 of the strange quark 
mass, although here Lg = 48 was required. Thus, we see domain wall fermions producing 
the desired light surface states with small mixing, even for relatively strong couplings. 

We have measured hadron masses and for lattice scales 1 GeV < <2 GeV and have 
studied scaling in this region. Our determinations of /,r involve not only mres as mentioned 
above but also the measurement of the Z-factor for the local axial current. We find f-j^/mp 
evaluated at the ruf + m^es = point to be scaling very well, while for mjq /nip the scaling 
violations may be at the 6% level. However, scaling seems at least as good as that seen for 
staggered fermions at similar lattice spacings and similar to that found for Wilson fermions 
with a clover term |Q . This is in accord with general expectations that finite lattice spacing 
errors will enter domain wall fermion amplitudes at O(a^) p!9| , |88 |. 



41 



Our results demonstrate that quenched domain wall fermions do exhibit the desired good 
chiral properties, even at relatively strong couplings. The residual quark mass effects, which 
break the full global symmetries to leading order in a, can be eliminated by an appropriate 
choice of ruf, so that low energy physics should be well described by an effective theory 
with the continuum global symmetries. Quenched chiral logarithm effects may appear for 
quenched domain wall fermion simulations, as they do for other fermion formulations, but 
present no new difficulties. For large enough volumes, the effects of topological near-zero 
modes are suppressed and the small rrif region can be investigated. For larger values of rrif, 
where these zero mode effects are suppressed by the quark mass, one has a formulation of 
lattice QCD with the full global symmetries realized to order and an effective quark mass 
of rrif + rrij-es- Thus, the domain wall formulation provides a powerful new tool which can be 
used, even within the quenched approximation, to study many of the outstanding problems 
in particle and nuclear physics for which chiral symmetry plays an important role. 

Note added: After this paper was essentially complete, the recent work of the CP- 
PACS collaboration became available The reader is referred to this paper for another 



discussion of some of the topics presented here. 
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TABLES 



TABLE L Simulation parameters for the quenched results obtained using the hybrid Monte 



Carlo method. The mass ranges referred to are specified in Table IIL The spectrum column 
contains the number of configurations on which hadron mass measurements were performed while 
the {qq) column shows the number of configurations used to compute the chiral condensate. Finally, 
within parenthesis in the last column we specify the number of random noise sources (hits) that 
were used in each of these {qq) measurements. All of the calculations described in this table used 
the domain wall height parameter M5 = 1.65. 
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TABLE II. Simulation parameters for the quenched results obtained using the heatbath 
method. The column labeled sweeps records the number of Monte Carlo sweeps between suc- 
cessive measurements. The remaining notation is the same as that used in Table |. 
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u 
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0. 
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o X 61 


A 

64 


l.b5 


r nnn 
oUOU 
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U.Uz— U.ZZ 


/b 


(J 


0. lU 


1 ^3 QO 

ib X oz 


O A 

Z4 


l.bo 


c^nnn 
oUOU 


n no n oo 
U.Uz— U.ZZ 


^Q 

/ o 


70(1) 


0. lU 


1 fi3 QO 

ib X 61 


QO 
OZ 


i.bO 


onnn 
zUUU 


n nnnoc^ n nno n nn n n/i 
U.UUUzO— U.UUo, U.UU— U.U4 


n 
U 




0. lU 
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ib X 61 


4o 


i.bO 
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U.UZ— U.ZZ 
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U 

^!1 
Di 
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0. 
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40 
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iUb 
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U 


0.00 
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on 
2U 


i.y 


c^nnn 
OUUU 


n no^ n nvc; 
U.UzO— U.U ( 


1 nn 
iUU 
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U 


O.oO 


q3 w qo 
X oZ 


QO 
OZ 


i.bO 




n nni nni n no n i n 
U.UUi— U.Ui, U.Uz— U.iU 


n 
U 


onn 
zUU 


O.oO 


1 f;3 V "^9 
iD X oZ 


QO 
OZ 


i.bO 


1 nnn 

iUUU 


n nni nni n no n i n 

U.UUi— U.Ui, U.UZ— U.iU 


n 

U 


01 ^^ ^ 


D.U 


io 


ib 


i.O 


onnn 
ZUUU 




QO 
OZ 




D.U 


1 fi3 QO 

ib X OZ 


1 o 

iz 


1 Q 
i.O 


onnn 
ZUUU 


n no 
U.Uz 


Ob 


n 
U 


D.U 


1 ^^3 w QO 

ib X oz 


ib 


i.O 


p^nnn 
OUUU 


n ni n n/i 
U.Ui— U.U4 


oO 


85(1) 


D.U 


1 fi3 w QO 

ib X oZ 


iD 


i.O 


onnn 
ZUUU 


n nnn 
U.UUU 


oi fi 

Zib 


n 
U 


D.U 


1 ^^3 w QO 

ib X oZ 


iD 


i.O 


onnn 
ZUUU 


n nni 
U.UUi 


ooo 
zzy 


n 
U 


6.0 


16^ X 32 


16 


1.8 


2000 


0.02 


56 





6.0 


16^ X 32 


16 


1.8 


2000 


0.00025-0.008, 0.00-0.04 





120(1) 


6.0 


16^ X 64t 


16 


1.8 


2000 


0.01-0.05, 0.075, 0.1, 0.125 


98 





6.0 


16^ X 32 


24 


1.8 


2000 


0.01-0.04 


76 





6.0 


16^ X 32 


24 


1.8 


2000 


0.02 


56 





6.0 


163 X 32 


24 


1.8 


2000 


0.00025-0.008 0.00-0.04 





110(1) 


6.0 


16^ X 32 


32 


1.8 


2000 


0.02 


72 





6.0 


16^ X 32 


48 


1.8 


2000 


0.02 


64 






t This extent of 64 in the time direction was achieved by "doubling" Nt = 32 lattice con- 
figurations in the time direction so the resulting gauge field background has an unphysical 
t — > t + 32 periodicity. 
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TABLE III. Here we list the explicit masses that are included in the various mass ranges 
referred to in Tables | and ||. 



mass range 


mass values 


0.00-0.04 


0.00, 0.005, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035 and 0.04 


0.00025-0.008 


0.00025, 0.0005, 0.001, 0.002, 0.004 and 0.008 


0.001-0.01 


0.001, 0.004, 0.007 and 0.01 


0.01-0.04 


0.01, 0.015, 0.02, 0.025, 0.03, 0.035 and 0.04 


0.01-0.05 


0.01, 0.02, 0.03, 0.04, 0.05 


0.02-0.20 


0.02, 0.04, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.18 and 0.20 


0.02-0.22 


0.02, 0.06, 0.10, 0.14, 0.18 and 0.22 


0.02-0.10 


0.02, 0.04, 0.06, 0.08 and 0.10 


0.025-0.075 


0.025, 0.0325, 0.05, 0.0625 and 0.075 



TABLE IV. Results for hadron masses at (3 = 5.7, 8^ x 32, M5 = 1.65, = 10 from 87 
configurations. 



ruf 




mp 


rriN 


O.Oz 


O.o28(lzj 


0.82(3) 


1.24(6) 


0.04 


0.604(12] 


0.874(19) 


1.32(4) 


0.06 


0.076(14) 


0.922(18) 


1.40(3) 


0.08 


0.744(12) 


0.971(15) 


1.48(3) 


0.10 


0.808 10 


1.017(12) 


1.56(2 


0.12 


0.869(9) 


1.063(10) 


1.63(2) 


0.14 


0.929(8) 


1.108(9) 


1.705(18) 


0.16 


0.987(5) 


1.153(8) 


1.779(15) 


0.18 


1.043(5) 


1.199(7) 


1.853(14) 


0.20 


1.097(5) 


1.244(6) 


1.928(13) 


TABLE V. 


Results for hadron masses at /3 = 


5.7, 8^ X 32, M5 = 1.65, 


= 16 from 67 


configurations. 








m/ 




rup 


ruN 


0.02 


0.483(18) 


0.87(12) 


1.25(14) 


0.04 


0.562(12) 


0.89(6) 


1.36(10) 


0.06 


0.635(7) 


0.93(4) 


1.44(8) 


0.08 


0.702(5) 


0.98(2) 


1.52(7) 


0.10 


0.768(5) 


1.019(18) 


1.59(6) 


0.12 


0.831(5) 


1.064(15) 


1.65(4) 


0.14 


0.892(6) 


1.109(13) 


1.71(4) 


0.16 


0.954(5) 


1.153(12) 


1.78(3) 


0.18 


1.012(5) 


1.199(10) 


1.85(3) 


0.20 


1.072(6) 


1.243(10) 


1.93(2) 
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TABLE VI. Results for hadron masses at /? = 5.7, 8^ x 32, Mg = 1.65, = 24 from 84 
configurations. 



mf nip uin 

002 0.44(2) 0.84(8) 1.4(4) 

0.06 0.613(12) 0.94(3) 1.31(6) 

0.10 0.756(8) 1.016(15) 1.50(3) 

0.14 0.882(5) 1.102(11) 1.662(16) 

0.18 0.999(4) 1.189(9) 1.817(12) 

0.22 1.108(4) 1.278(8) 1.959(11) 



TABLE VII. Results for hadron masses at /3 = 5.7, 8^ x 32, M5 = 1.65, Ls = 32. The 

results for m/ = 0.02 arc obtained from 278 configurations, those for m/ = 0.04 are from 184 
configurations, while the others are from 94 configurations. 

ruf m-jr TUp rriN 

002 0.405(6) 0.83(5) 1.17(11) 

0.04 0.502(5) 0.87(4) 1.16(5) 

0.06 0.595(9) 0.92(2) 1.36(8) 

0.10 0.743(7) 0.995(16) 1.50(3) 

0.14 0.872(6) 1.082(11) 1.66(2) 

0.18 0.991(5) 1.178(8) 1.822(16) 

0.22 1.104(4) 1.274(7) 1.970(15) 



TABLE VIII. Results for hadron masses at /3 = 5.7, 8^ x 32, M5 = 1.65, Ls = 48. The results 
for ruf = 0.04 are obtained from 335 configurations, while the others are from 169 configurations. 
One of the original 336 configurations had eigenvalues very close to zero requiring nearly 11,000 
conjugate gradient iterations to converge. The resulting pion propagator was so large as to dom- 
inate the average of the (7r"(a;)7r"(0)) correlator for the mf = 0.0 case. We omitted this single 
coiiiigurHiiou from (liis Mruilysis. 



rUf rrtTT rup ruN 

002 0.374(10) 0.99(12) 1.07(19) 

0.04 0.490(4) 0.87(2) 1.22(6) 

0.06 0.580(7) 0.95(3) 1.38(5) 

0.10 0.730(5) 1.016(13) 1.52(2) 

0.14 0.860(4) 1.098(8) 1.658(17) 

0.18 0.981(4) 1.184(6) 1.809(14) 

0.22 1.093(4) 1.272(5) 1.962(12) 
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TABLE IX. Results for hadron masses at /? = 5.7, 8^ x 32, Mg = 1.65, = 64 from 76 
configurations. 







nip 


VHN 


0.02 0.364(14) 




0.98(17) 


1.2(3) 


0.06 0.563(8) 




0.96(5) 


1.18(4) 


0.10 0.719(8) 




1.01(2) 


1.42(3) 


0.14 U.854(7j 




l.U89(12j 


l.bZyZ) 


0.18 0.978(6) 




1.176(9) 


1.784(17) 


0.22 1.097(5) 




1.265(7) 


1.938(14) 


A., rtesuits lor naaron masses 


1 at p 


— 0.1, ID X oZ, IVl^ 


— i.DO, ijg — Zft irom to 


configurations . 








nif rriTT 




rup 


rriN 


0.02 0.412(5) 




0.85(4) 


1.17(9) 


0.06 0.597(4) 




0.909(12) 


1.29(4) 


0.10 0.743(4) 




0.990(9) 


1.47(2) 


0.14 0.873(4) 




1.082(6) 


1.642(14) 


0.18 0.994(3) 




1.175(5) 


1.800(16) 


0.22 1.105(4) 




1.265(4) 


1.953(16) 


IA-JdIjIij A.1. liesults tor nadron masses 


at P = 


= 5.7, 16^ X 32, Ms = 


1.65, Lg = 48. The results 


lor T/if — u.uz, U.U41:, u.uD, u.uo ana u.iu are 


obtained from 106 configurations, while the others are 


from 6-1- (■< )iiIl^Q,'urMi;k)iis. 








rrif rriTr 




rrip 


niN 


0.02 0.383(4) 




0.88(5) 


1.00(4) 


n n/i n /I so/zi ^ 

U.U4 U.4oZ(^4J 




0.86(2) 


1.23(3) 


0.06 0.577(2) 




0.918(18) 


1.260(15) 


0.08 0.650(3) 




0.951(9) 


1.375(16) 


0.10 0.729(2) 




0.989(9) 


1.449(10) 


0.14 0.865(3) 




1.083(7) 


1.623(14) 


0.18 0.986(3) 




1.173(5) 


1.786(12) 


0.22 1.097(4) 




1.264(5) 


1.942(12) 
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TABLE XII. Results for hadron masses at ^ = 5.85, 12^ x 32, M5 = 1.9, Lg = 20, 100 config. 



mf 






nip 


mN 


0.0250 


0.359(4) 




0.627(15) 


0.844(24) 


0.0375 


0.426(3) 




0.650(10) 


0.919(16) 


0.0500 


0.483(3) 




0.675(8) 


0.985(13) 


0.0625 


0.536(3) 




0.706(6) 


1.044(10) 


0.0750 


0.585(3) 




0.738(5) 


1.10(3(9) 


TABLE XIII. 


Results for hadron masses at /3 = 


6.0, 16^ X 32, M5 = 1.8 


, Lg = 16 from 85 


configurations. 
















nip 


nij\f 


0.010 


0.203(3) 




0.442(10) 


0.621(30) 


0.015 


0.239(3) 




0.451(7) 


0.648(21) 


0.020 


0.270(3) 




0.462(6) 


0.668(15) 


025 


0.298(3) 




0.475(5) 


0.686(12) 


0.030 


0.324(2) 




0.488(5) 


0.706(10) 


0.035 


0.348(2) 




0.502(4) 


0.729(9) 


0.040 


0.371(2) 




0.515(4) 


0.752(9) 


TABLE XIV. 


Results for hadron 


masses with nondegenerate valence quarks at P = 6.0, 


16^ X 32, M5 = l.i 


3, Ls = 16 from 98 configurations. 






mf[l) 


m/(2) 






nip 


0.010 


0.020 




0.238(2) 


0.441(12) 


0.030 


0.020 




0.298(2) 


0.471(8) 


0.040 


0.020 




0.325(2) 


0.487(7) 


0.050 


0.020 




0.349(2) 


0.502(7) 


0.075 


0.020 




0.406(2) 


0.538(6) 


0.100 


0.020 




0.457(2) 


0.574(5) 


0.125 


0.020 




0.505(2) 


0.608(4) 
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TABLE XV. Results for hadron masses at /3 = 6.0, 16^ x 32, M5 = 1.8, = 24 from 76 
configurations. 



nif 




nip 




0.010 


0.201(5) 


0.423(11) 


0.664(43) 


0.015 


0.236(4) 


0.441(9) 


0.644(25) 


0.020 


0.267(3) 


0.458(7) 


0.653(18) 


0.025 


0.295(3) 


0.473(6) 


0.674(15) 


0.030 


0.320(3) 


0.487(5) 


0.698(12) 


0.035 


0.345(3) 


0.502(5) 


0.723(11) 


0.040 


0.368(2) 


0.516(4) 


0.748(10) 



TABLE XVL 


Valence extrapolations (a + bmf) for ml, m 


9 and mN at P = 5.7, 


8^ X 32, 


M5 = 1.65. Tlie fitting ranj 


;es used are described in Section [II C 






mass 


Ls 




b 


/dof 


2 

mi 


10 


0.201(17) 


4.55(17) 


5.35 


rup 


10 


0.792(19) 


2.26(6) 


0.55 


rriN 


10 


1.18(4) 


3.74(17) 


0.13 


2 


16 


0.147(32) 


4.32(36) 


1.59 


nip 


16 


0.798(24) 


2.23(10) 


0.17 


TUN 


16 


1.21(8) 


3.60(29) 


0.28 


ml 


24 


0.093(16) 


4.82(12) 


1.47 


rrip 


24 


0.786(21) 


2.24(7) 


0.34 




24 


1.17(3) 


3.60(14) 


1.01 


mi 


32 


0.076(10) 


4.83(12) 


4.81 


rUp 


32 


0.753(26) 


2.36(11) 


0.43 


ruN 


32 


1.14(5) 


3.78(22) 


1.14 


ml 


48 


0.042(8) 


4.90(9) 


0.012 


rrip 


48 


0.790(13) 


2.18(5) 


0.50 


rriN 


48 


1.13(4) 


3.79(15) 


0.73 


ml 


64 


0.039(11) 


4.75(21) 


6.28 


rUp 


64 


0.756(22) 


2.28(8) 


0.65 


mN 


64 


1.03(6) 


4.06(22) 


3.23 
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TABLE XVII. Valence extrapolations (a + brnj) for m^, rrip and m^v at /3 = 5.7, 16^ x 32, 
M5 = 1.65. The fitting ranges used are described in Section [II C. 



mass 


Ls 


a 






b 


xVdof 


ml 


24 


0.072(6) 






4.83(9) 


7.26 


nip 


24 


0.764(12) 






2.27(5) 


1.10 


ruN 


24 


1.10(3) 






3.87(15) 


0.69 


nip 


48 


0.775(18) 






2 2Q(7) 


1 04 


rriN 


48 


1.03(4) 






4.13(17) 


3.58 


TABLE XVIII. Valence extrapolations (a + bmf) for 


ml, 


mp and m^ at (3 


= 5.85, 12^ X 32, 


M5 = 1.9. The 


masses from all five values of are included 


in the fits. 




mass 


Ls 


a 






b 


xVdof 


ml 


20 


0.024(3) 






4.27(4) 


2.8 


rrip 


20 


0.549(14) 






2.50(17) 


n ss 

U.oo 


rriN 


20 


0.74(2) 






4.9(3) 


n SI 

U.Oi 


TABLE XIX. Valence extrapolations (a + bmf) for 


ml, 


mp and m^ at (3 


= 6.0, 16^ X 32, 


M5 = 1.8. The fitting ran 


ges used are described in Section 


inc. 




mass 


Ls 


a 






b 


xVdof 


2 

m^ 


16 


0.0098(20) 






3.14(9) 


0.029 


rrip 


16 


0.404(8) 






2.78(11) 


0.48 


ruN 


16 


0.566(21) 






4.66(29) 


0.34 


ml 


24 


0.0094(26) 






3.09(7) 


0.32 


rrip 


24 


0.400(10) 






2.86(12) 


0.38 


rriN 


24 


0.546(19) 






5.05(33) 


0.65 


TABLE XX. Results for fits of {qq) to the form given in Eq. |4^. 


X NtX Ls 


P 


a_i ao 






ai 




8^ X 32 X 32 


5.7 


6.0(6) X 10-^ 1.76(3) X 10~ 


-3 




6.53(4) X 10-2 


4.0(4) X lO^^ 


8^ X 32 X 48 


5.7 


6.8(7) X IQ-^ 1.92(5) x 10" 


-3 




6.04(14) X 10^2 


1.7(2) X 10-^ 


16^ X 32 X 32 


5.7 


2.5(4) X 10"^ 1.86(2) x 10" 


-3 




6.53(2) X 10-2 


9.3(9) X 10-3 


16^ X 32 X 16 


6.0 


1.0(1) X IQ-^ 3.87(8) X 10" 


-4 




8.64(1) X 10^2 


5.6(3) X lO"'' 


16^ X 32 X 24 


6.0 


9.1(10) X lO"'^ 3.62(9) X 10" 


-4 




8.64(2) X 10^2 


1.1(1) X 10""^ 
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TABLE XXI. Results for m,r from different correlators for ruf < 0.01. There should be a single 
value for the pion mass, determined at asymptotically large times, irrespective of the correlator 
used. Since the correlators give different masses for the fitting ranges used, the localized topo- 
logical near-zero mode effects are important. Here PP, AA and PP-I-SS represent the correlators 
(^''(x)7r'*(0)), (^g(x)^g(O)) and (7r'^(x)^'*(0)) + {a{x)a{0))c respectively. 



p 


V 


Ls 


correlator 


mj 






xVdof 


5.7 


8^ X 32 


48 


PP 


0.0 




0.273(39) 


1.3 ±0.6 


5.7 


8^ X 32 


48 


AA 


0.0 




0.197(11) 


1.5 ± 1.0 


1^ 7 
0. / 


O A 


^0 


-T r -TOO 


u.u 




U. IZol^Zi J 


1.0 It U.o 


5.7 


163 X 32 


48 


PP 


0.0 




0.200(5) 


1.9 ±0.9 


5.7 


16^ X 32 


48 


AA 


0.0 




0.193(7) 


1.6 ±0.9 


5.7 


16^ X 32 


48 


PP-FSS 


0.0 




0.191(7) 


1.1 ±0.6 


6.0 


16^ X 32 


16 


PP 


0.0 




0.151(8) 


0.6 ±0.6 


6.0 


16^ X 32 


16 


AA 


0.0 




0.098(7) 


0.6 ±0.5 


6.0 


16^ X 32 


16 


PP+SS 


0.0 




0.017(43) 


1.0 ±0.6 


6.0 


16^ X 32 


16 


PP 


0.001 




0.141(6) 


0.9 ±0.6 


6.0 


16^ X 32 


16 


AA 


0.001 




0.118(6) 


1.1 ±0.9 


6.0 


16^ X 32 


16 


PP+SS 


0.001 




0.082(11) 


0.5 ±0.4 


TABLE XXIL Results for 


residual mass at /3 = 5.8. 


5 and 6.0. The /? = 5.85 calculation was 


performed 


on a 12^ X 32 with anti-periodic boundary conditions, M5 = 


1.9 and the ratio R{t) from 


Eq. 77 averaged over the time range 6 < t < 26. The /3 


= 6.0 calculation, described 


in the text, 


was performed on a 16^ 


X 32 lattice with M5 = 1.8. 












Ls 




nif 


# of config 






5.85 


20 




0.05 




100 




0.00281(8) 


6.0 


12 




0.02 




56 




0.00239(6) 


6.0 


16 




0.02 




56 




0.00124(5) 


6.0 


24 




0.02 




56 




0.00059(4) 


6.0 


32 




0.02 




72 




0.00044(4) 


6.0 


48 




0.02 




64 




0.00027(3) 
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TABLE XXIII. Results for residual masses at /3 = 5.7, M5 = 1.65. 



Lattice size 


Ls 


nif 


^ of config. 




8^ X 32 


32 


0.02 


184 


0.0106(2) 


8^ X 32 


32 


0.04 


184 


0.0105(2) 


8^ X 32 


48 


0.04 


335 


0.00688(13) 


16^ X 32 


48 


0.02 


50 


0.0072(9) 


16^ X 32 


48 


0.04 


50 


0.0071(4) 


16^ X 32 


48 


0.06 


50 


0.0066(6) 


16^ X 32 


48 


0.08 


50 


0.0065(5) 


16^ X 32 


48 


0.10 


50 


0.0063(4) 



TABLE XXIV. Results for fits to the form predicted for a quenched chiral logarithm, Eq. 82 
for the 16'^ x 32 simulations at /? = 5.7 with Lg = 48. For comparison, we have also included 
the column Xjntcpt which gives the x intercepts predicted by the simple linear fits of Eq. 75 in 



the heavier mass range 0.02 < nif < 0.08. Here PP, AA and PP+SS represent the correlators 
(7r"(x)7r'*(0)), (ylg(x)ylg(O)) and (7r"(x)7r"(0)) + {a{x)a{0))c respectively. 



Correlator 



ao 



ai 



a2 



xVdof 



S^intcpt 



PP 

AA 
PP+SS 



4.7(3) 
4.1(3) 
4.3(3) 



0.0085(7) 
0.0073(10) 
0.0075(11) 



-0.008(28) 
-0.07(4) 
-0.05(4) 



1.9 ±4.7 
3.6 ±2.4 
4.8 ±3.2 



-0.0092(12) 
-0.0108(7) 
-0.0104(7) 



TABLE XXV. Results for Za at f3 = 6.0, 16^ x 32, M5 = 1.8. 



Lg rrif configurations Za 

12 002 56 0.7560(3) 

16 0.02 56 0.7555(3) 

24 0.02 56 0.7542(3) 

32 0.02 72 0.7535(3) 

48 0.02 64 0.7533(3) 
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TABLE XXVI. Results for the correlator amplitude A and at /? = 5.7, 8^ x 32, M5 = 1.65, 
Lg = 32, 94 configurations. Parameters with subscript AA are obtained from the axial vector 
current correlator. Parameters with subscript PP are obtained from the pseudoscalar density 
correlator. Za = 0.7732 and rrires = 0.0105 are used in the /^r calculation as described in the text. 



rrif Aaa iU)AA ^PP {U)pp 

0.02 0.0077(11) 158(9) 0.176(18) 145(8) 

0.06 0.0132(13) 172(8) 0.145(10) 173(6) 

0.10 0.0197(15) 188(7) 0.143(7) 194(4) 

0.14 0.0267(15) 202(5) 0.151(6) 214(4) 

0.18 0.0346(15) 216(4) 0.164(5) 233(3) 

0.22 0.0432(16) 229(4) 0.180(5) 252(3) 



TABLE XXVIL Results for the correlator amplitude A and at /3 = 5.7, 8^ x 32, M5 = 1.65, 
Ls = 48, 169 configurations. Parameters with subscript AA are obtained from the axial vector 
current correlator. Parameters with subscript PP are obtained from the pseudoscalar density 
correlator. Za = 0.7732 and mres = 0.00688 are used in the /jr calculation as described in the text. 



rrif Aaa {U)aa ^pp {U)pp 

0.02 0.0058(5) 134(5) 0.202(13) 143(7) 

0.06 0.0114(7) 152(4) 0.144(7) 162(4) 

0.10 0.0172(8) 167(3) 0.138(5) 180(3) 

0.14 0.0240(9) 182(3) 0.146(4) 198(3) 

0.18 0.0317(10) 196(3) 0.159(4) 216(2) 

0.22 0.0403(12) 209(3) 0.176(5) 234(3) 



TABLE XXVIII. Results for the correlator amplitude A and at /? = 5.7, 16^ x 32, M5 = 1.65, 
Lg = 48. Parameters with subscript AA are obtained from the axial vector current correlator. 
Parameters with subscript PP are obtained from the pseudoscalar density correlator. Za = 0.7732 



and iriy 


= O.OOGcScS i\rv used iu (lie , 


[j^ calcuLiiiou as d(^sc 


riV)('d iu the text. 




ruf 


Aaa 


(Maa 


App 


{f-^)pp 


0.00 


0.0028(2) 


133(5) 


0.34(4) 


136(11) 


0.02 


0.0063(3) 


142(3) 


0.180(10) 


138(4) 


0.04 


0.0092(5) 


153(4) 


0.140(6) 


150(4) 


0.06 


0.0115(4) 


156(3) 


0.139(6) 


163(3) 


0.08 


0.0147(5) 


167(3) 


0.127(4) 


169(3) 


0.10 


0.0174(5) 


171(2) 


0.136(5) 


182(3) 


0.14 


0.0248(10) 


188(4) 


0.147(8) 


201(5) 


0.18 


0.0326(12) 


202(4) 


0.161(7) 


220(4) 


0.22 


0.0411(15) 


214(4) 


0.178(7) 


239(4) 
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TABLE XXIX. Results for the correlator amplitude A and at /? = 6.0, 16^ x 32, M5 = 1.8, 
Lg = 16, 85 configurations. Parameters with subscript AA are obtained from the axial vector 
current correlator. Parameters with subscript PP are obtained from the pseudoscalar density 
correlator. Za = 0.7555 and mres = 0.00124 are used in the /jr calculation as described in the text. 



rrif Aaa iU)AA ^PP {U)pp 

0.010 0.00100(13) 144(10) 0.050(6) 149(9) 

0.015 0.00132(14) 153(8) 0.040(4) 151(8) 

0.020 0.00161(15) 159(8) 0.036(4) 156(8) 

0.025 0.00189(16) 164(7) 0.034(3) 161(8) 

0.030 0.00216(17) 168(7) 0.032(3) 166(8) 

0.035 0.00243(18) 172(6) 0.032(3) 171(8) 

0.040 0.00269(18) 175(6) 0.031(3) 176(7) 



TABLE XXX. Results for the correlator amplitude A and at /3 = 6.0, 16^ x 32, M5 = 1.8, 
Lg = 24, 76 configurations. Parameters with subscript AA are obtained from the axial vector 
current correlator. Parameters with subscript PP are obtained from the pseudoscalar density 
correlator. Za = 0.7542 and mres = 0.00059 are used in the /jr calculation as described in the text. 



ruf 


Aaa 


{U)aa 


App 


(/7r)pp 


0.010 


0.00089(14) 


138(11) 


0.043(6) 


134(10) 


0.015 


0.00110(11) 


141(7) 


0.034(4) 


138(8) 


0.020 


0.00130(11) 


144(6) 


0.031(3) 


144(6) 


0.025 


0.00152(11) 


148(5) 


0.029(2) 


149(6) 


0.030 


0.00176(12) 


153(5) 


0.0280(18) 


154(5) 


0.035 


0.00201(14) 


157(5) 


0.0275(16) 


159(5) 


0.040 


0.00227(15) 


162(5) 


0.0272(15) 


164(5) 



TABLE XXXI. Linear fit parameters, /t^, fx and fK/fn determined from the axial vector 
current correlator and the pseudoscalar density correlator. 



p 


V 


Ls 


correlator 


intercept 


slope 


U 


fK 


Ik/U 


5.7 


8^ X 32 


32 


axial 


152(9) 


352(38) 


148(10) 


162(8) 


1.094(16) 


5.7 


8^ X 32 


32 


pseudoscalar 


142(6) 


505(25) 


137(6) 


157(5) 


1.146(14) 


5.7 


8^ X 32 


48 


axial 


130(4) 


364(18) 


127(4) 


145(4) 


1.142(11) 


5.7 


8^ X 32 


48 


pseudoscalar 


135(4) 


453(21) 


132(4) 


154(4) 


1.171(13) 


5.7 


16^ X 32 


48 


axial 


136(3) 


362(33) 


133(4) 


149(2) 


1.122(14) 


5.7 


16^ X 32 


48 


pseudoscalar 


129(3) 


525(36) 


125(4) 


149(2) 


1.188(18) 


6.0 


16^ X 32 


16 


axial 


138(10) 


958(195) 


137(11) 


156(8) 


1.134(37) 


6.0 


16^ X 32 


16 


pseudoscalar 


138(10) 


938(235) 


137(10) 


155(8) 


1.131(40) 


6.0 


16^ X 32 


24 


axial 


128(10) 


847(222) 


127(10) 


143(7) 


1.124(41) 


6.0 


16^ X 32 


24 


pseudoscalar 


123(10) 


1031(232) 


122(11) 


142(7) 


1.156(49) 
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TABLE XXXII. ResTilts for a ^ using nip extrapolated to + mj-es = to set the scale and 
for mN/nip extrapolated to this same value of m/. 



_^ V Ls g-^GeV) ruN/mp 

"5?7 8^ X 32 32 1.058(40) 1.51(5) 

5.7 8^ X 32 48 0.994(18) 1.45(3) 

5.7 16^ X 32 48 1.013(24) 1.40(5) 

5.85 12^ X 32 20 1.419(34) 1.34(5) 

6.0 16^ X 32 16 1.922(40) 1.42(4) 

6.0 16^ X 32 24 1.933(50) 1.38(4) 



TABLE XXXIII. A variety of expressions for {qq) in lattice and physical units. The quantity b 
in the final column comes from our earlier (a+bnif) fits to m^. The final column gives a convention- 
ally normalized value of the chiral condensate which is to be compared with the phenomenological 

value of (229 ±9)^ 



X NtX Ls 




oo 


{QQ)mf =—mies 


&/^/48 




8^ X 32 X 32 


5.7 


1.76(3) 10-3 


1.07(3) 10-3 


1.7(2) 10-3 




8^ X 32 X 48 


5.7 


1.92(5) 10-3 


1.50(5) 10-3 


1.8(1)10-3 




16^ X 32 X 16 


6.0 


3.87(8) 10-^ 


2.80(9) lO-'' 


3.3(5) 10-4 


245(7)MeV 


16^ X 32 X 24 


6.0 


3.62(9) 10-4 


3.11(10) 10-4 


2.6(5) 10-4 


256(8)MeV 
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FIGURES 
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FIG. 1. (qq) for quenched simulations done on 8^ x 32 lattices (o) and 16^ x 32 lattices (□) 
at f3 = 5.7 with Lg = 32. The smaller volume shows a pronounced rise as m/ ^ as is expected 
if unsuppressed zero modes are present. For the larger volume, the effect of topological near zero 
modes is reduced if not eliminated. This is expected since {\iy\)/V should fall as l/W. 
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FIG. 2. The pion mass squared versus nif from (7r"(x)7r°(0)) (□), {Aq{x)Aq{0)) (o) and 
(7r''(x)7r"(0)) + (cj(x)cr(0))c (O) for quenched simulations done on 8^ x 32 lattices at /3 = 5.7 
with Lg = 48. For mj = 0.0, the correlators all give different masses due to the differing topologi- 
cal near-zero mode contributions for each one. For large enough x, all the correlators should give 
the same mass. However, this limit requires a large volume which is expected to suppress such 
zero- mode effects. The dotted line is the fit of Eq. p9, the solid line is from Eq. 7C and the dashed 



line is from Eq. 71 
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FIG. 3. The pion mass squared versus rrif from (7r'*(x)7r"(0)) (□), {Aq[x)Aq{^)) (i 
(7r"(x)7r'^(0)) + ((T(x)cr(0))c (O) for quenched simulatfons done on 16^ x 32 lattices at /3 

The star is the value of measured from Eq. |77| and its error bar in the 



and 
5.7 



with 



48. 



horizontal axis is too small to show on this scale. The solid line is the fit to the (^Q(a;)^Q(0)) 
correlator for = 0.02 to 0.08 given in Eq. 75, while the dotted line is for the (7r'*(x)7r'*(0)) 
correlator for = 0.0 to 0.08 as given in Eq. 7C. 
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FIG. 4. The dependence of the residual mass for 16^ x 32 lattices at /3 = 6.0. The long-dashed 



line is the fit given in Eq. [7q, the short-dashed line is the fit from Eq. 79 and the solid line is the fit 



given in Eq. 81. Each of the three fits is made to all of the points shown. We have employed an 
intermediate non-pGrturbcttive rcnorinctlization to convert the plotted values of tti^qs into the MS 
scheme at // = 2 GeV. 
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FIG. 5. Results for at /? = 6.0 with a 16^ x 32 lattice and Lg = 16 plotted as a function of 

mj. The open circles are obtained from the {Aq{x) Aq{{^)) correlator, while the open diamonds are 
obtained from the (7r''(x)7r"(0)) correlator. We also show the linear fits which are used to determine 
our estimate for /jr and Jk- The vertical dashed lines identify the values for m/ which locate the 
chiral limit, m/ = — mres and give the physical ratio for rriK/mp. The solid symbols represent the 
extrapolations to the point ruf = — rrires and interpolations to the kaon mass. 
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FIG. 6. Effective mass, m^sit) is plotted for the 16^ x 32, /3 = 6.0, = 16, m/ = 0.01 
calculation of the vr, p and nuclcon masses. While plateau regions for t > 5 arc easily identified for 
uItj- and nip, the nucleon fit is less satisfactory. Although a plateau may be recognized for 5 < i < 8, 
the rapidly growing errors make such an identification problematic for this case. More satisfactory 
nucleon plateaus are seen for the larger values of m/. 
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FIG. 7. Distribution of and the mass difference {mi{t) — mi{tanu))/cri for t > tmin, evaluated 
for three vahies of tmin for the case of the p mass and a 16'^ x 32 lattice, with /3 = 6.0 and TUf > 0.01 
and Ls = 12, 16, 24, 32 and 48. Both distributions appear reasonable for each value of imin with 
only small improvement as tmin increases from 5 to 9. We choose to quote values of nip for all these 
cases using the value tmin = 7. 
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FIG. 8. The p (o) and nucleon masses (□) plotted as a function of for the case of /? = 5.7, 
M5 = 1.65, Ls = 32 and a 8^ x 32 lattice. The lines represent least squares fits whose parameters 



appear in Table XVI while the data plotted appears in Table VII. Note the relatively low value 
for the nif = 0.04 nucleon point results from the comparison of two somewhat different data sets. 
As mentioned in the text, the linear fit was obtain from the 94 configurations identified in Table | 
while the nif = 0.02 and 0.04 points plotted also include the further 184 configurations referenced 
in Table 0. 
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FIG. 9. The p (o) and nucleoli masses (□) plotted as a function of for the case of /3 = 5.85, 
M5 = 1.9, Lg = 32 and a 8^ x 32 lattice. The lines represent least squares fits whose parameters 



appear in Table XVIII while the data plotted appears in Table XII 
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10. The p (o) and nucleon masses (□) plotted as a function of m/ for the case of /3 = 6.0, 
Ls = 16 and a 16'^ x 32 lattice. The lines represent least squares fits using the param- 



FIG. 
Ms = 1. 

eters appearing in Table pCIX| while the data plotted appears in Table [XIII| . In addition to these 
hadron masses computed for the case of equal mass quarks, we have also plotted the p mass for 
the case of non-degenerate quarks given in Table pCIV as a function of the average quark mass, 
(mj(l) + mj(2))/2. These points are plotted as filled diamonds. 
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FIG. 11. (qq) for quenched simulations done on 8^ x 32 lattices at /? = 5.7 for 



0.05 



32 



and Lg = 48 (□). The more pronounced rise as m/ ^ for Lg = 48 shows that the expected 
topological near-zero modes have smaller values for Aj and/or 5mi for this larger Lg. 
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FIG. 12. {qq) for quenched simulations done on 16^ X 32 lattices at /3 = 6.0 for Ls = 16 (o) 
and Lg = 24 (□). The more pronounced rise as ruf ^ ioi Lg = 24 shows that the expected 
topological near-zero modes have smaller values for and/or Srrii for this larger Lg. 
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FIG. 13. Evolutions of —{qq) (solid line) and — (9759) (dotted line) for 16^ x 32 lattices at 
f3 = 6.0 with Lg = 16. For smaller values of mf the evolutions show pronounced fluctuations which 
have opposite sign for (qq) and {q^5q), indicating the presence of eigenfunctions of Dh at this Lg 
which are very good approximations to the exact topological zero modes expected at Lg 
Note that the vertical scale increases for the smaller values of m/. 
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FIG. 14. The quantity m^(m/ = 0) from (7r"(x)7r"(0)) for quenched simulations done on 8^ x 32 
lattices at /3 = 5.7 versus Lg. This graph is an updated version of an earlier result based on 
part of the data presented here. While a slow decrease in m^imf = 0) as Lg increases from 32 
to 64 is now visible, the effect is much less dramatic than the drop seen in the more accurate 
values of mres, which decrease from 0.0105(2) to 0.0071(4) as Lg increases from 32 to 48. This 
contrast presumably results from the effects of both zero- modes and non-linearity for m^(m/) in 
the quenched approximation. 
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FIG. 15. The pion effective mass as a function of the source-sink separation, t, for 8^ x 32 
lattices at /3 = 5.7 with Lg = 48 and m/ = 0.0. The upper panel is from {Aq{x)Aq{0)) (o), the 
middle from (7r"(x)7r"(0)) (□) and the lower from (7r"(a;)7r"(0)) + {a{x)a{0))c {<>)■ The nif = 0.0 
points in Figure come from fitting from t = 7 to t = 16. 



73 



0.55 

0.5 

0.45 
0.55 

0.5 

0.45 
0.55 

0.5 

0.45 







OAA 



o o 



□ PP 



[] 



O PP+SS 



<> 



^ 5 X ^ 



8 
t 



12 



16 



FIG. 16. The pion effective mass as a function of the source-sink separation, t, for 8^ x 32 

lattices at /3 = 5.7 with Lg = 48 and nif = 0.04. The upper panel is from {Aq{x)Aq{0)) (o), the 
middle from (7r"(x)7r"(0)) (□) and the lower from (7r"(x)7r"(0)) + {a{x)a{0))c (O). The effective 
masses from the different correlators are quite consistent. 
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FIG. 17. The evolution of point source correlators at t = 8 for 8^ x 32 lattices at f3 = 5.7 with 

Lg = 48 and rrif = 0.0 The upper panel is (7r"(x)7r"(0)), the middle —{a{x)a{0))c and the lower 
{Aq{x)Aq{0)) . The large fluctuations that are common to (-7r"(a:)7r"(0)) and —{a{x)a{0))c are due 
to zero modes and show that they dominate the ensemble average for the correlator at this t. 
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FIG. 18. The evolution of point source correlators at t = 8 for 8'^ x 32 lattices at /3 = 5.7 with 
Ls = 48 and m/ = 0.04 The upper panel is (7r"(x)7r"(0)), the middle —{a{x)a{0))c and the lower 
{Aq{x)Aq{0)). Zero mode effects seem entirely absent from these evolution plots. 
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FIG. 19. The pion mass squared versus irif from (7r'*(x)7r"(0)) (□), {Aq{x)Aq{0)) (o) and 
(7r"(a;)7r"(0)) + {a{x)a{0))c {O) for quenched simulations done on 16'^ x 32 lattices at (3 = 6.0 
with Ls = 16. For ruf = 0.0 and 0.001, the correlators all give different masses due to the differing 
topological near-zero mode contributions for each one. For larger m/, the pion mass determination 
from (7r"(x)7r"(0)) + ((7(x)(t(0))c is likely contaminated by the heavy mass states in the {a{x)a{0))c- 
The dotted line is the fit of Eq. [f3, the solid line is from Eq. 73 and the dashed line is from Eq. 74. 
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FIG. 20. The pion effective mass as a function of the source-sink separation, t, for 16^ x 32 
lattices at /3 = 6.0 with = 16 and nif = 0.001. The upper panel is from {Aq{x)Aq{0)) (o), the 
middle from (7r"(a;)7r"(0)) (□) and the lower from (7r"(x)7r"(0)) + {a{x)a{0))c (O). 
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FIG. 21. The pion effective mass as a function of the source-sink separation, t, for 16^ x 32 

lattices at /3 = 6.0 with Lg = 16 and mj = 0.01. The upper panel is from {Aq{x)Aq{0)) (o), the 
middle from (7r"(x)7r"(0)) (□) and the lower from (7r"(x)7r"(0)) + {a{x)a{0))c {O). 
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FIG. 22. For each mj, the average value of rri-Tr is calculated for the three correlators and the 
graph above shows the deviation of each correlator from the average. For each nif, the result from 
{A^ix)A^{Q)) is shifted slightly to the right and the result from (7r"(x)7r"(0)) + (c7(x)<t(0))c to the 
left for clarity. No systematic deviation is visible from the data. 
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FIG. 23. The pion effective mass as a function of the source-sink separation, t, for 16^ x 32 
lattices at /3 = 5.7 with Lg = 48 and nif = 0.0. The upper panel is from {Aq{x)Aq{0)) (o), 
the middle from (7r"(a;)7r"(0)) (□) and the lower from (7r"(x)7r"(0)) + {a{x)a{0))c (O). All three 
correlators give reasonable effective masses and the fitted masses agree. 
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FIG. 24. The evolution of point source correlators aX t = 8 for 16^ x 32 lattices at /? = 5.7 with 
Lg = 48 and mj = 0.0 The upper panel is (7r"(a;)7r"(0)), the middle —{a{x)a{0))c and the lower 
{Aq{x)Aq{0)) (AA). There are few, if any, contributions from topological near-zero modes. 
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FIG. 25. The per degree of freedom for linear fits of versus mj for 16^ x 32 lattices at 
P = 5.7 with Ls = 48 from {A^{x)A^{0)) (o), (7r''(x)7r"(0)) (□) and (7r"(x)7r"(0)) + {a{x)a{0))c 
{O). Only the range ruf = 0.02 to 0.08 gives a fit with an acceptable value for per degree of 
freedom. After using the large volume to eliminate zero modes, and presumably also finite volume 
effects, we have evidence for a non-linear dependence of on m/. 



83 



0.015 



0.01 



CO 
0) 



0.005 







[] 



[] 



n [] 







T_ n 








L : 


] 



[] 



0.0105(2) 



0.00688(13) 







10 
t 



15 



20 



FIG. 26. The residual mass for 8^ x 32 lattices with Lg = 32,48 and m/ = 0.04 at /3 = 5.7. 
The labels for the horizontal lines are the averages over the range 4 < t < 16 with jackknife errors. 
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FIG. 27. The pion mass squared versus m/ from (7r"(a;)7r''(0)) (□), {Aq{x)Aq{0)) (c 
(7r''(x)7r'*(0)) + {a{x)a{0))c i^) for quenched simulations done on 16^ x 32 lattices at f3 
with Lg = 48. The star is the value of mres as measured from Eq. ^ and its error bar in the hori- 
zontal axis is too small to show on this scale. The solid line is a fit of the {Aq{x)Aq(0)) correlator 
for nif = 0.0 to 0.08 to the quenched chiral logarithm form given in Eq. This fit gives the pion 
mass vanishing in very good agreement with the value of rwres determined from Eq. 77. 
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FIG. 28. A three dimensional "Lego" plot showing the matrix elements of between all 
nineteen eigenvectors found for one of the better configurations in our sample of 32, evaluated at 
mf = 0. The height of the box located by horizontal coordinates represents the magnitude 
of the matrix element (Aj:f^j|F^| A//j). The five zero modes, all nearly eigenvectors of F^ with 
eigenvalue +1, are easily identified. The remaining seven pairs are also very evident corresponding 
to the expected \Ah) and F^|A//) = | — Ah) eigenstates. 
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FIG. 29. The distribution of the 4-dimensional norm, M{s) for the first two zero modes shown 
in Figure ^ as a function of s. Note both states are tightly bound to the s = Lg — 1 waU, as are 
the other three zero modes states. 
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FIG. 30. The distribution of the 4-dimensional norm, M{s) for the first pair of non-zero modes 



shown in Figure 28 as a function of s. 
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FIG. 31. The distribution of values of the quadratic fit parameter 6mi defined in Eq. 48. 
These parameters were determined from a total of 576 eigenvalues obtained from 32 configurations 
computed with /? = 6.0, 16^ and Lg = 16. The peak of the distribution lies remarkably close to 
the value we find for the residual mass, rrij-es = 0.00124(5). 
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FIG. 32. The renormalization constant Za obtained for a 16^ x 32 lattice with ^5=16 and 
/3 = 6.0 (o) and that for a 8^ x 32 lattice with Ls=48 and ^ = 5.7 (O). The labels for the 
horizontal lines are the averages, with jackknife errors, over the ranges 4 < t < 14 and 18 < t < 28. 
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FIG. 33. Results for /^r at /? = 5.7 with a 8^ x 32 lattice and Lg = 48 plotted as a function oiruf. 
The open circles are obtained from the axial vector current correlator, while the open diamonds are 
obtained from the pseudoscalar density correlator. Wc also show the linear fits which are used to 
determine our estimate for fj^ and /k- The vertical dashed lines identify the values for rrif which 
locate the chiral limit, nif = — mres and give the physical ratio for niK/mp. The solid symbols 
represent the extrapolations to the point ruf = —mres and interpolations to the kaon mass. 
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FIG. 34. Results for at /? = 5.7 with a 16^ x 32 lattice and Lg = 48 plotted as a function of 

nif. The open circles arc obtained from the axial vector current correlator, while the open diamonds 
are obtained from the pscudoscalar density correlator. We also show the linear fits which are used 
to determine our estimate for and fx- The fits are done to the points with nif = 0.02 — 0.10. 
The vertical dashed lines identify the values for mf which locate the chiral limit, ruf = —rrires and 
give the physical ratio for mK/mp. The solid symbols represent the extrapolations to the point 
ruf = —rrires ^^nd interpolations to the kaon mass. 
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